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abstract 

The technique of differential very-long-baseline 
interferometry (VKBX) has been used to ' 12( 

^r e irirana° f 17 t Xnar Tand™tes with uncertainties 
ilia tAan 0"005 of geocentric arc. These measurements 
mi vielded improved determinations of the selenodetrc 
coordinates^o^the Apollo landing sites, and of the phy- 
sical libration of the moon. 

RV means of a new device, the Differential Doppler 
Receiver (DDR) . instrumental errors were -ducedjo.les^ 

of“the NASA°Spacef light Tracking ^m^lonT”* 
i qttvm ) and used in an extensive program of observati 
(STDN ) ana. __ na( . a oT-rt-ained over a 16— month 

y u ^ i— 

“rd-digSfhLmlnifcoplicients of thl 

SS-SSSM 1 ^- VLB! data and from the cons rstency 

radia^coordinate^and STS rse coordinates. 

Values determined for the libration parameters 
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8 G f and C.,., have uncertainties smaller than the uncertain- 
ties 3 Stained when laser ranging data alone is used in the 
solution. The rms of the postfit residuals for the VLBI 
observations is 16° of phase (at 2,3 GHz), about 2 times 
larger than the random noise level. The systematic compo- 
nents in the residuals may result from unmodeled propagation- 
medium effects. 
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CHAPTER I 
Introduction 

The moon has always been a prime object of study 
by astronomers, and observations of its motion have im- 
portant uses in many fields, for example to test gravita- 
tional theories and to measure variations in the earth's 
rotation. In order to interpret these observations, 
theories of both the orbital and the librational motions 
of the moon are required. However, until recently the 
development of libration theories [ 20 , 23] received much 
less attention than did orbital theories, due to the 
limitations of ground-based passive optical observations, 
which for a long time were the only observations available. 
The departure from uniform rotation represented by the 
forced or "physical" libration is less than 2 seconds of 
arc as viewed from earth, near the resolution limit im- 
posed by "seeing" fluctuations. Recently, however, the 
placement on the lunar surface of man-made optical re- 
flectors and radio transmitters and the development of 
instrumentation allowing extremely accurate earth-based 
measurements have stimulated a new interest in the theory 
of lunar librations. This thesis is concerned with the 
first use of differential very-long-baseline interferometry 
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(VLBI) to measure the relative positions of radio trans- 
mitters on the lunar surface, and to determine the phy- 
sical libration by observation of the apparent time- 
variation of the relative transmitter positions. 

Radio tracking of objects on or near the moon has 
been carried out in the U.S. since the first Ranger mis- 
sion in 1964 and has continued through the Surveyor, 

Lunar Orbiter, and Apollo programs. However, the track- 
ing transmitters for these missions operated only for 
short times, and the tracking data were not suitable 
for the determination of positions on the moon with uncer- 
tainties less than 100 meters [44, 36, 29] . The po- 

tential for a longer-duration and higher-accuracy program 
of radio observations was created by the last five Apollo 
missions, each of which left on the moon a nuclear-powered 
ALSEP* with an S-band (~2.3 GHz) transmitter (Figure 1.1) • 
These transmitters have only free-running crystal oscil- 
lators and are not able to be tracked by the usual two- 
way Doppler-shift or ranging methods used for spacecraft 
with transponders. However, they do provide strong sig- 
nals suitable for angular-position measurements by VLBI. 


*ALSEP is an acronym for Apollo Lunar Surface Experiments 
Package . 
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Fig. 1.1 ALSEP transmitters operating in 
May 1975. 

In interferometry, the phases of the signals re- 
ceived from a source simultaneously at two separated 
ground stations are compared. As the interferometer 
baseline rotates relative to the source, the difference 
between the phases changes. The direction of the source 
with respect to the baseline is determined from this phase 
variation (see Figure 1.2). The existence of more than 
one ALSEP transmitter permits the use of the technique of 
differential interferometry [11], in which the difference 
in phase of the signal received from one ALSEP at two 
stations is subtracted from the corresponding dif- 
ference for another ALSEP. The resulting doubly- 
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differenced observable is sensitive to the relative 
right ascension and declination of the two ALSEP trans- 
mitters, and is relatively insensitive to those other 
parameters such as the moon's orbital position, the 
observing site positions, receiver local oscillator stabi 
lity and the phase fluctuations introduced by the earth's 
atmosphere and ionosphere, which tend to affect the obser 
vations of both transmitters equally. 



Fig. 1.2 Radio source position determination by VnBI. 

A change of phase by 2 tt rad (1 full cycle) cor 
responds to a source position angle change of 
A/D rad, where A is the wavelength of the radio 
signal and D is the distance between the stations 
projected onto a plane perpendicular to t e 
direction of the source. For ALSEP observa- 
tions such as reported in this thesis 
(A *13.2 cm, D * 4000 km) the phase change 
can be determined with an uncertainty of less 
than 10° . Were there no systematic errors present, 
this phase uncertainty would correspond to an angular 
uncertainty of less than 10” 9 rad or a displacement un- 
certainty of only 40 cm at the distance of the moon 119] . 








-14- 


The idea to apply the VLBI technique to study the 
motions of the moon stems from a 1968 proposal made to 
the National Aeronautics and Space Administration (NASA) 
by H. F. Hinteregger and I. I. Shapiro of through 

the M.I.T. Center for Space Research, to add a wideband 
(- 50 MHz) noise source, center frequency near 8 GHz, to 
each ALSEP package. The choice of such an X-band fre- 
quency range was based primarily on the desire for the 
reduction of ionospheric effects on the VLBI observables 
to a negligible level and for the achievement of compati- 
bility with then-existing ground-based VLBI equipment. 

The scientific objectives included the determination of 
the motion of the moon's center of mass through differen- 
tial VLBI observations of the ALSEPs and compact extra- 
galacti - radio sources that are nearly occulted by the 
moon, ls well as the determination of the moon's motion 
about its center of mass through differential VLBI ob- 
servations of two or more ALSEPs. This Lunar Radio Beacon 
proposal was not accepted by NASA. Thereafter, C. C. Counsel 
man III, of M.I.T., along with Hinteregger and Shapiro, 
elicited the cooperation of I. M. Salzberg of the Goddard 
Space Flight Center of NASA to attempt to use the narrow- 
band, lower frequency, S-band signals from the ALSEP tele- 
metry system to meet many of the same scientific objectives. 


-15- 


The resultant MIT-GSFC cooperation led to NASA's Space- 
flight Tracking and Data Network (STDN) undertaking 
differential VLBI observations of the ALSEPs in early 
1971, as soon as two ALSEP transmitters were operating 
on the moon. The accuracy of these observations was 
degraded by instrumental errors equivalent to positional 
uncertainties on the lunar surface of the order of 100 m, 
due mainly to the fact that the signals from the two 
ALSEPs did not pass through identically the same receiving 
equipment and were sampled separately before a difference 
observable was formed (see Figure 1.3) . 

In order to reduce these instrumental errors, Hin- 
teregger and Counselman designed and built a new device, 
called a Differential Doppler Receiver (DDR) . Used with 
the antennas and S-band receivers of the STDN, the DDR 
extracts the carrier signals of two ALSEPs from the 
intermediate-frequency passband of a single receiver and 
mixes them together. The resulting difference frequency, 
after a constant "bias" is subtracted, is multiplied by 
a factor of 360 to improve phase resolution, then counted 
digitally. A differential VLBI observable is obtained 
by numerical subtraction of the counts recorded simul- 
taneously at two stations. 
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Fig. 1.3 An example of the results of early dif- . 

ferential interferometric observations 
of ALSEP telemetry transmitters. Note 
the presence of phase noise on the order 
of a cycle peak to peak and drifts of 
several cycles, corresponding to tens of 
meters displacement on the lunar surface. 


In October 1972, Hinteregger and I installed pro- 
totype models of the DDR in the STDN tracking stations 
at Merritt Island, Florida, and Goldstone, California, 
and conducted successful observations of the ALSEP 12 
and ALSEP 14 transmitters during a six-hour tracking 
period. Analysis of the data from this first experiment 
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indicated that measurement errors had been reduced to the 
equivalent of less than a meter of displacement on the 
moon [13] (Figure 1.4). 


POST-FIT RESIDUALS 28 OCT I9T2 
ALSEPS 12 AND 14 

MERRITT ISLAND FLA,- GQLDSTONE, CAL IE 



Fiq. 1.4 Residuals (observed minus computed values) 

' for differential interferometric phase 

observations of ALSEP 12 and ALSEP 14, ex- 
pressed directly as fractions of a cycle at 
S-band (wavelength, 13.2 cm) and as the equiva 
lent (projected) displacement of an ALSEP 
transmitter on the surface of the moon. 


In early 1973 we installed two DDRs at each of six 
stations of the STDN (Figure 1.5). Since that time, 
simultaneous observations have been made of two pairs of 

ALSEPs an average of twice per week on each of two 
earth baselines. The data set analyzed in this thesis 


consists of the observations during the 
16-month period between March 1973 and July 1974 



Fig. 1.5 


Stations of the NASA Spacecraft Tracking and 
Data Network that have engaged m track 

inq of the ALSEPs are Ascension Island (ACN) ; 
Madrid , Spain (MAD); Merritt Island, Florida 
(MIL); Corpus Christ!, Texas (TEX); Goldstone, 
California (CDS) ; and the Network Test and 
Training Facility at Greenbelt, Maryland (ETC) 


Analysis of these data has been carried out using 
the M.I.T. Planetary Ephemeris Program (PEP), a large 
computer program originally developed at the M.I.T. 

under the direction of Michael E. Ash 


Lincoln Laboratory 
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and Irwin I. Shapiro. PEP is able to utilize data from 
many types of astronomical observations simultaneously 
to estimate various physical parameters [2,6 ] . 

Recently, Martin A. Slade, in his M.I.T. thesis research, in- 
corporated into PEP an accurate, numerically-integrated 
model for the lunar orbit [36]. The use of this model 
has been essential to my work. I have made the further 
modifications to PEP necessary to process counted-cycle 
differential VLBI observations, and to model more ac- 
curately the moon's physical libra tion. 

Before and during the period of our observations, 
members of the NASA Lunar Laser Ranging Experiment (LURE) 
team [7 ] have collected and analyzed laser ranging ob- 
servations from the McDonald Observatory in Texas to 
the retroref lectors on the moon (principally the re- 
flectors at the Apollo 11, 14, and 15 landing sites). 
Since differential VLBI observations have little direct 
sensitivity to the orbital motion of the moon, I have 
used laser data to improve PEP's lunar ephemeris. In 
addition, I am indebted to members of the LURE team for 
the development of much-improved models of the lunar 
physical libration. In my data reduction I have used 


- 20 - 


principally a numerically integrated model developed 
at the Jet Propulsion Laboratory by Williams et al . 

[40] . The JPL numerical integration derives its initial 
conditions by fitting to a semi-analytic model developed 
by D. H. Eckhardt of the Air Force Cambridge Research 
Laboratories [ 15,16 ] . 

The ALSEP relative coordinates determined from these 
VLBI observations should be useful in the definition of 
a unified selenodetic coordinate system that incorporates 
both the high accuracy of laser ranging and the wide cover- 
age of Apollo orbital photography and earth-based tele- 
scopic photography. In particular, the VLBI observations 
serve to extend to the Apollo 12, 16, and 17 landing 
sites the positional accuracy obtained from the laser 
ranging observations of the Apollo 11, 14, and 15 sites. 
The Apollo 16 and 17 coordinates obtained by VLBI are 
the most important because these landing sites are 
covered by the metric photographs taken by the orbiting 
command module during these two missions. The Apollo 15 
laser ranging retroref lector has been located relative 
to the ALSEP transmitter at that sit© with an uncertainty 
of less than one meter by L. A. Schimerman et al. [30] at the 
Defense Mapping Agency Aerospace Center, using photographs 
taken on the lunar surface by the astronauts. 


- 21 - 


In the following chapters I try to document in some 
detail each major aspect of the ALSEP VLBI experiment. 

I begin with a basic discussion of counted-cycle VLBI 
and of the sensitivity of such observations to relative 
ALSEP positions. In Chapter III, a description is given 
of the Differential Doppler Receiver and its use with 
the NASA STDN Unified S-Band tracking system. Chapter IV 
deals with the software developed to analyze the data, 
including the special processing algorithms required for 
the doubly-differenced observable, the documentation of 
the model used for the lunar rotation, and the method of 
parameter estimation. The chapter concludes with an 
analysis of the major sources of error which might de- 
grade the results. The results of my solutions for the 
ALSEP coordinates and the libration parameters are given 
in Chapter V. A concluding chapter suggests the likely 
gains from further analysis of these data and from the 
application of the differential VLBI technique to other 
problems in dynamical astronomy. 
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CHAPTER II 
Basic Concepts 

2 . 1 Counted-cycle VLBI 

The usual VLBI system consists of two or more stations 
each with four key elements: 1) a high-gain antenna 

with a low-noise radio frequency (RF) amplifier; 2) one 
or more stages of mixing with "local oscillators" to 
convert the received RF signals to low "video" frequencies 

3) some means of recording the converted signals; and 

12 

4) a stable (1 part in 10 or better) frequency standard 
to provide a reference for the local oscillators and a 
time base for the recording. Systems differ in the ways 
in which the interferometric phase is determined from 
the recordings made at the two stations . In the con- 
ventional technique, the signals are recorded directly 

on magnetic tape in either analog or digital form at 
each station. The tapes are later collected, played 
back together, and the data cross-correlated to determine 
the differential phase delay. For observations of natural 
radio sources, which emit essentially random signals, 
statistical cross-correlation of the signals is necessary 
to determine the phase difference. This technique also 
offers the possibility of determining the group delay 
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difference with useful accuracy if the bandwidth of the system 
is sufficient. 

For observations of artificial radio transmitters which 
emit phase-coherent signals, an alternative to the relatively 
expensive process of direct signal recording and cross-correla- 
tion is available. An oscillator can be locked in phase to the 
signal received at each station and a digital counter used con- 
tinuously to count the cycles of this oscillator. At suitable 
intervals the numerical value of the continuously-accumulating 
count can be recorded*. The oscillator can be locked to a 
multiple, n, of the signal frequency, so that each count 
represents 1/nth cycle of the signal phase. Thus, a phase 
resolution of a degree or less is possible. An interferometric 
phase observable is formed simply by subtracting the counts 
corresponding to the same time at the two stations. For 


4r 

Counted-cycle VLBI observations of spacecraft can easily be 
made between any pair of sites equipped to perform Doppler 
tracking. Ondrasik and Rourke of JPL have analyzed the tech- 
nique as a supplement to the two-way Doppler tracking of inter- 
planetary spacecraft [25]. Although some observations were 
made of Mariner 9 [ 45] , application of this technique has been 
limited by the lack of sufficiently stable frequency standards 
at most sites: Atomic hydrogen maser oscillators are required 

for useful results. If two spacecraft are observable simul- 
taneously, the differential VLBI technique can be employed and 
the commonly used cesium-beam standards would be adequate 
since the effects of local-oscillator frequency variations would 
tend to cancel. 


differential observations involving signals from two trans- 
mitters, the oscillator may be locked to the difference between 
the frequencies of the two signals received at a given station. 

In present practice, the phase delay observable determined 
by counting cycles differs from that determined by direct sig- 
nal recording in two ways. In direct recording VLB I the re- 
ceived signals are compared for common transmission times, since 
in the cross-correlation one of the recordings is shifteid with 
respect to the other by an amount equal to the propagation time 
difference, determined from prior knowledge. In all of the cycle- 
counting systems employed thus far, the counts have been sampled 
simultaneously at the two receiving sites, so that the phase dif- 
ference determined by subtraction refers to signals transmitted 
at different times. One consequence of this transmission-time 
difference is that short-term (time scale of milliseconds) phase 
fluctuations of the transmitter that are not filtered out by the 
phase-locked loops may add noise to the measurement, which can- 
cels in the interferometric observable only when the propagation 
times to the two stations are equal. An additional problem with 
sampling simultaneously at both receivers is the possible intro- 
duction of systematic errors into the model of the observable due 
to persistent (over minutes or hours) drifts in the transmitter 
frequency. This effect in the ALSEP observations will be analyzed 
in Section 4.5. In principle, both of these errors could 
be reduced by offsetting the receiver 
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sampling times during observations using a priori in- 
formation about the differential delay between stations. 

The systematic effect of persistent transmitter frequency 
drift may also be accounted for in the theoretical formu- 
lation of the counted-cycle observable, as described in 
Section 4.2. 

The second difference between the two VLBI methods, as 
now implemented, involves the technique of conversion of the measured 
phase difference to a corresponding delay difference. This conversion, 
which is necessary in order to relate the observed quan- 
tity to the geometry of the observing and transmitting sites, involves 
the use of the observed frequency. In conventional VLBI ob- 

servations there are two possible cases. For a continuum 
source, there is no characteristic source frequency, and 
the observed frequency is determined by the character- 
istics of the receiver: usually, although not necessarily, 

by the total local oscillator (LO) frequency. For a 
monochromatic source, however, there is a characteristic 
frequency associated with the received signal. This 
signal is recorded within the video band at whatever 
frequency the source appears relative to the total LO. 

In this case, the conversion of the measured phase to a 
delay should be based on a suitably defined signal frequency. 
If, instead, the measured phase were referred to the LO fre- 
quency using the group delay, the derivative of phase with 
respect to frequency, then a phase 
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error would be introduced, equal to the product of the 
error in the value of group delay used, and the offset 
of the LO frequency from the signal frequency. In a 
cycle-counting VLBI system, no extrapolation to refer 
the phase from the received frequency to the LO frequency 
at a station is necessary, since the local oscillator 

is variable and continuously tracks the phase, and there- 
fore the frequency, of the received signal. In practice 
we convert phase to delay using an estimate of the 
actual transmitted frequency, obtained from a measurement 
of the received frequency at one site and a priori know- 
ledge of the Doppler shift of the signal received at 
that site. This procedure is discussed in detail in 


Section 4.2. 


Compared with direct-recording VLBI, cycle counting offers 
significant simplifications in observation procedure and data 
processing. Numerical data are obtained directly in real time 
and can be recorded at a relatively low sampling rate, for 
example, once per minute. (In fact, the data could be written 
down by hand, using pencil and paperl) In addition, the dif- 
ferential phase observable may be formed in real time by com- 
municating via the telephone or teletype the numbers recorded 
at each station, thus allowing the experimenter to check that 
valid data are being obtained. 


rr sesra* 




The principles underlying the use of differential 
VLB! observations to determine the relative coordinates 
of a pair of ALSEP transmitters can be understood by means 
of a simplified geometric analysis in which the moon is 
assumed to be infinitely far from the earth and to remain 
stationary during a single day's tracking period*. The 
geometry of a pair of stations observing a single such 
ALSEP is shown in Figure 2.1 where £>{t) is the baseline 
vector and is the unit vector in the direction of 
ALSEP i. The difference between the phases 



Figure 2.1 Geometry for two stations observing 


a source at infinity. 


Ir 

The approximations introduced in this section are not made in the actual 
data processing. See Chapter IV. 
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of the signals received simultaneously at the two stations, 
usually called the interferometric "fringe phase", is 
simply 

4>i (t) = £ b(t) ■d i (2.2.1) 

where w is the angular frequency of the transmitted 
signal, c is the speed of light, and is the path- 

length difference. The two vectors may be written 

S(t) = B{cos D [cos A(t)l + sin A(t)5] + sin D k} 

( 2 . 2 . 2 ) 

= cos 6^ (cos a i 1 + sin a ± 5) + sin k (2.2,3) 

where B is the baseline length; i, 5, R, are unit vectors, 

1 in the direction of the vernal equinox, & in the direction 
of, the earth’s angular velocity vector, and j orthogonal 
to i and k so as to complete a right-handed system; A, 

D, and a^, 6^ refer to the right ascension and the decli- 
nation of the baseline and the source, respectively. 

If two ALSEP transmitters are observed simultaneously, 
the differential fringe phase is given by 

Acj> (t) = 4> 1 (t) - 4> 2 <t) = | S(t) * (d x - S 2 ) (2.2.4) 

where I have assumed that the two transmitter frequencies 
are the same. (All of the ALSEPs do in fact have differ- 
ent frequencies, but the frequency differences are small 
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and do not need to be considered here.) In order to 
express A<J> (t) in terms of the small differential right 
ascension, Aa, and differential declination, AS, between 
the directions to the two ALSEPs, we write 

a 2 = a l + 

62 = + A 6 

and substitute these expressions in Equations (2.2.3) and 

(2.2.4). Since the actual geocentric angular separation 

-3 

of the ALSEPs is less than 5 x 10 radians, we expand 
the trigonometric functions in powers of Aa and A 6 , ig- 
noring terms of order higher than the first in these angles. 
Expanding yields 

^1 “ ^2 - [AS sinS^ cosa-^ + Aa cos S-^ sina^]£ 

+ [ASsinS-^ sina-^ - AacosS^ cosa^JJ 

- AS cosS 1 k (2.2.5) 

Substituting Equations (2.2.2) and (2.2.5) in 
Equation (2.2.4) gives 

A<f>(t) = ™{AS [cosDsinScos (a - A(t)) - sinDsinS] 

+ Aa [cosDcosSsin (a - A(t))]} (2.2.6) 

in which the subscript on a and S has been dropped. 

By counting cycles of the DDR outputs and taking 
the difference between the counts at the two stations we 


-31- 


can measure the change in the differential fringe phase 
during the period of observation, but not the "absolute" 
value at a given time. Hence the constant term in 
Equation (2.2.6) is not observable. The observable quantity 


Ad) . (t) - — cosD [ A5 sinScos(a - A(t) ) 

Y obs c 

+ AacosSsin(a - A(t))] + A 4> 


(2.2.7) 


where Ad> is an unknown constant which accounts both for 
o 

the constant term in Equation (2.2.6) and for the 
difference between the arbitrary initial values of the 
cycle counters at the two stations. 

Equation (2.2.7) shows that the coefficients of the 
sinusoidal and cosinusoidal variations of Ac ^ obs ^ are 
proportional, respectively, to Aa and AS. Thus, determining 
Aa and A 6 for an ALSEP pair corresponds to determining the 
amplitudes of the quadrature (sine and cosine) components 
of the variation of Acj> obs (t) . If the observation period 
is long enough to allow separate determination of both 
quadrature components, or, equivalently, of both the 
amplitude and the phase of the diurnal sinusoid, then both 
Aa and AS may be estimated with the uncertainties given 


a Aa ” a s wBcosDcosS 


J c mBcosDsinS 


( 2 . 2 . 8 ) 

(2.2.9) 


where cr and a are the uncertainties in the deter- 
s c 

mination of the sine and cosine coefficients, respectively, in 
Equation (2.2.7). Typical values with present equipment 
for the parameters in Equations (2.2.8) and (2.2.9) 

9 

are a ~ a c ~ 20° , w ~ 2ir x 2.3 x 10 radians/sec and 
BcosD (the equatorial projection of the baseline) 

~ 3500 km. If the moon is at a relatively high decli- 
nation, say 6 = 20°, then ~ 0"0005 and 

°A<5 ~ 0V0012. Of course, at low lunar declinations, 

A6 is poorly determined. 

These values provide only crude estimates of the 
uncertainties of the determinations of Aa and A6. In 
practice , due to the problem of separating the sine and 
cosine coefficients from the constant component in 
Equation (2.2.7), the formal standard errors obtained 
for a particular observation period depend strongly on the 
time span of the observations and on the characteristics of 
the earth baseline. To investigate these problems, 

Frankston [18 ] has performed a series of three-parameter 
least-squares solutions based on the model of Equation (2.2.7) 
using simulated observations for each of the available STDN 
baselines (see Figure 1.5). For each solution he assumed 
that the error in each observation was 10 ° and that ob- 
servations were made at one- hour intervals throughout 
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the "observation window" for a particular baseline. This 
window was defined as the time during which the moon was 
at least 20° above the horizon for both stations. For 
each observation period three parameters were estimated: 
Aa, AS, and the constant, A4> Q - For all baselines which 
had equatorial projections longer than 1500 km and ob- 
servation windows longer than 3 hours, the solutions 
yielded formal standard errors for Aa between 070003 and 
0700 6 , equivalent to 0.5 m to 10 m in distance at the 
moon. For AS, considering only observations made when the 
declination of the moon was at least 10° (about three- 
fourths of each month) , he found that the formal standard 
errors ranged from 070004 to 07014 or 0.8 m to 30 m 


in equivalent distance on the surface of the moon. 

The type of baseline most useful for determining 
each of the relevant parameters can be deduced by in- 
spection of Equation (2.2.7). Since observations typically 
are obtained over only 4-6 hours, the curvature of the 
diurnal-period sinusoid is rather poorly determined from 
the observed change in differential fringe phase. The 
best-determined quantity is the average slope of A<j)(t); 
consequently , the quadrature component of A<j> (t) which is 
best determined is the one which has its maximum slope 
during the period of the observations. The cosine function 
multiplies A 6 and has its maximum slope when 
a - A(t)= ±90°. This implies that for the best deter- 
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mination of M, the baseline should be perpendicular 
to the direction of the moon at the midtime of the obser- 
vations, a condition best satisfied by a baseline that 
is oriented east-west on the earth. Conversely, the sine 
function which multiplies Aa has its maximum slope when 
a - A(t) = 0 or 180% so that, for the best determination 
of Aa, the baseline vector should be parallel to the 
meridian plane of the moon at the midtime of the observa- 
tions. This condition is satisfied by a baseline that is formed 

by a pair of stations with nearly equal longitudes, but that 
has a large equatorial projection. 

From the foregoing analysis, it is clear that ob- 
servations should be made with both east-west and north- 
south baselines in order to obtain optimum determination of 
both Aa and A6 from the relatively short observing periods 
available. The application of this strategy in the con- 
text of station availability and other scheduling con- 
straints is discussed in Section 3.4. 

By combining observations made on different days, 
when the moon presents different aspects to the earth, 
we can use this differential VLBI technique to determine 
the third (radial) component of the separation, or baseline 
vector, between the ALSEPs. This determination is possible 
because the moon’s geometric, or "optical” libration, which 
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varies by ±7° in both longitude and latitude, allows the 
interferometer to view the lunar baseline from different per- 
spectives during the course of a month. The standard error 
associated with the radial component is greater by a fac- 
tor of about 3 than the formal errors associated with the 
transverse components. 

The continuation of VLBI observations over a sufficient- 
ly long time also enables the moon ' s physical libration to 
be determined, through the variations of Act and A <5 which 
the libration produces. The apparent change in the 
relative angular position of a pair of ALSEPs produced by 
the libration depends on the orientation and length of 
the ALSEP baseline and on the nature of the libration . 

The libration may be characterized by rotations about three 
axes, one along the line of sight from earth to moon and 
two perpendicular to this line. Differential VLBI ob- 
servations are most sensitive to the librations about the 
line of sight, since the ALSEPs are all located relatively 
close to a single plane normal to the line of sight. 

ALSEPs 12 and 17, for example, are separated 
by apprQximately 1300 km along a nearly 
east-west line (see Figure 1.1). Determination of AS 
with an uncertainty of O'.'OOl for this pair corresponds to 
a determination of the libration about the line of sight 
with an uncertainty of 0'.'001 x ~ i 3 qq ^ " 0 1*3 , where 
384000 km is the earth-moon distance. The sensitivity 
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of VLBI observations to libration about the axes per 
pendicular to the line of sight is proportional to the 
component of an inter-ALSEP baseline along the line of sight. 

The baseline which has the greatest component in this 
direction is formed by ALSEP 14 (longitude -17°, latitude 
-4°) and ALSEP 17 (longitude 31°, latitude 20°)*. For 
this baseline, an uncertainty of 0V001 in the measurement of Aa 
is equivalent to an uncertainty in the libration in longi- 
tude (about a north-south axis) of 1 ' 1 6 , and an un- 
certainty in A 6 of 0 0 0 2 corresponds to an uncertainty of 
3V 2 in the libration in latitude (about an east-west axis) . 
Of course, the libration cannot be determined in all three 
degrees of freedom from one day's observations of a single 
pair of ALSEPs. Only two independent quantities, corres- 
ponding to Aa and AS, are observable. However, simul- 
taneous observations of a second ALSEP pair may be used to 
determine the third libration component. 

The assumption that the ALSEPs are located at an in- 
finite distance from the receiving stations has led to a 
relatively simple expression, Equation (2.2.7), which has 


helped us to see how a single day's observations can 


*The ALSEP 14 - ALSEP 17 pair is not directly observable 
with a single DDR since the frequency difference for this 
Dair is qreater than 3.3 MHz (see Section 3.2). However, 
the baseline may be synthesized in a solution by combining 
observations from two DDRs; for example one °DR observing 
ALSEP 14 with respect to ALSEP 12, the other ALSEP 17 with 

respect to ALSEP 12. 
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determine the angular separation of a pair of ALSEPs in 
both right ascension and declination. There is, however, 
additional information content in the observable provided 
by the lunar parallax. (The moon's parallax is defined 
loosely as the angle subtended at the distance of the moon 
by the earth's equatorial radius.) The parallax describes 
the extent to which an observer on the earth can view the 
moon from different perspectives during the course of a 
single day. For an interferometer we define the parallax 
as P/R, where R is the earth-moon distance and P is the 
distance of the baseline midpoint from the earth's axis 
of rotation. 

In Appendix A. 3 I derive an expression for A^^ft) for 
the two-dimensional case (D = 6 = 0) , including terms of 
first order in the parallax: 


A W (t) = A *o + i£ § £ ^{sin(a - A (t ) ) 

+ ^[sin 2 (a - A(t)) - j tan £ cos 2 (a - A(t))] 

+ Ji cos(aV 2 )cos 28 . (a . A(t)) 

r cos a 

- i-tan 2£ cos (a - A(t))]} (2.2.8) 

where r is the lunar radius, £ is the longitude of the 
m 

midpoint of the lunar baseline, and A£ is the separation of 
the ALSEPs in selenocentric longitude. Note that the first 
line of Equation (2.2.8) is exactly the same as 
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Equation' (2.2.7) for the case D = 6 = 0. 

The third and fourth lines in Equation (2.2.8) 
contribute no additional information since they 
have the same period as the main term*. However, the term 
in the second line, proportional to P/R, has a 12-hour 
period so that its signature is distinct from the signature 
of the main term. The determination of the amplitude and 
phase of the semi-diurnal sinusoid allows the determination 
of two additional parameters, P/R and £, for example. 
However, the uncertainties for these parameters will be 
larger than those for Aa by a factor of approximately 1/Aa 
which is about 400 for an ALSEP pair whose separation, 
when projected along the right-ascension direction, is 
about 1000 km. 


The simplified geometric analysis presented above is 
useful for the purpose of gaining qualitative insight into 
the way in which the observations determine the parameters 
of interest. However, it fails to answer quantitatively 
the question of what uncertainties may be expected if the 
various parameters are estimated from observations extend- 
ing over several months. The relatively complicated or- 
bital and rotational motions of the moon add significant infonration content 
to the data and decorrelate the estimates of parameters 
which could not be estimated 


The phas^ of the tem in the fourth line differs from the phase of the main 

[s™’Eq!StiOT r (2!2.7n : *** ““ Wl11 3130 haVe 3 °° S “ e 
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separately from a single day's observations. Clearly ad- 
vance knowledge of the parameter-estimation uncertainties 
obtainable from an experiment is important in planning an 
experiment. Therefore I used PEP to perform maximum- 
likelihood estimates of relevant parameters using simulated 
observations made on two days per week over the period 
March 1973 to July 1974 [19]. Each day's observations of 
an ALSEP pair consisted of only five differential- 
interferometric phase measurements, made one hour apart. 

The 1-a uncertainty of each measurement was assumed to be 
9°. I estimated simultaneously the six elements of the 
lunar orbit; the three selenocentric coordinates of each 
ALSEP; the three geocentric coordinates of each observing 
station; the two lunar moment-of-inertia ratios, 

3 [= (C-A)/B] and y [s (B-A)/C], upon which the charac- 
teristics of the physical libration mainly depend; and, 
for each daily tracking session, an additional parameter 
representing the arbitrary initial value of the differenced 
phase-counter readings. For the moon's orbital elements I 
assumed a priori 1-a uncertainties for a and e equivalent 
to 2.5 meters in radius, and for the angular elements, to 
100 m in the two orthogonal directions. The origin of the 
selenocentric coordinate system was defined by assuming 
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that the position of ALSEP 15 was known. The results of 
the simulation indicated that the coordinates of ALSEPs 12 
14/ 16, and 17 would be determined with 1-a uncertainties 
of approximately 1 meter in longitude and latitude and 3 
meters in radius. The uncertainty obtained for £ corres- 
ponded to an uncertainty of 0706 in 
the libration in latitude and 

node ; at one lunar radius, approximately the distance 
between the more widely separated ALSEPs , the corresponding 
displacement uncertainty is about 50 cm. The uncertainty 
obtained for y implies an uncertainty of 0723 in rotation 
about the spin axis (the libration in longitude) . The 
uncertainties obtained for the tracking station coordinates 
ranged from 5 to 15 meters in radius and longitude and 15 
to 20 meters in latitude. These uncertainties are com- 
parable to the uncertainties in the coordinates derived 
from other methods, so that the station coordinates should 
be included among the parameters to be estimated from the 
actual VLB I observations . 
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CHAPTER III 

Experimental Equipment and Procedure 

3.1 Unified S-Band System of the STDN 

The NASA Spaceflight Tracking and Data Network 
(STDN) was created in 1972 by the merger of the Manned 
Space Flight Network (MSFN) , which supported the Mer- 
cury, Gemini, and Apollo missions, and the Space Track- 
ing and Data Acquisition Network (STADAN) , which pro- 
vided primary tracking for the extensive NASA earth 
satellite program. Each of the former MSFN stations is 
equipped with either a 30-foot- or 85--foot-diameter 
parabolic antenna (or both) and an S -band tracking and 
telemetry receiving system capable of providing position 
determination and communications for spacecraft in the 
vicinity of the earth or moon. Since 1970 the network 
has used these stations regularly to observe the ALSEP 
transmitters on the moon, primarily for the purpose of 
demodulating and recording the information telemetered 
from the scientific experiments left on the surface by 
the astronauts. Because differential VLBI observations 
of the ALSEPs can be performed simultaneously with tele- 
metry support, an extensive program of VLBI observations 


has been conducted with little additional complication for 
the network . 

The central element of the Unified S-Band (USB) 
system is a Motorola superheterodyne receiver of a type 
built originally for use by the Deep Space Network* in the 
early 1960's. A simplified schematic diagram of the re- 
ceiver and its associated elements is shown in Figure 3.1. 
For differential Doppler tracking of ALSEPs, the S-band 
antenna is directed at the moon and the USB receiver is 
phase-locked to one ALSEP, which we call the "reference 
ALSEP". The carrier signal of the reference ALSEP appears 
in the center of the 10 ± 3.3 MHz intermediate 
frequency (IF) band of the USB receiver, at exactly the 
same frequency as, and with a fixed phase relationship to, 
the 10 MHz reference oscillator of the USB system. The 
signals from other ALSEPs will also appear within this IF 
band at various frequencies above and below 10 MHz, unless 
the radio frequency (RF) difference between the reference 
ALSEP and another exceeds about 3.3 MHz, in which case the 
latter signal falls outside the IF passband.** The USB 


*Operated for NASA by the Jet Propulsion Laboratory. 

'Jt'fc 

The antenna beams are sufficiently broad to receive signals from all 
ALSEPs simultaneously. 


S-BAND ANTENNA 


Reference ALSEP at f. 



Figure 3.1 


Schematic diagram of the STDN USB 
system. The signal of the reference 
ALSEP appears at 10 MHz in the out- 
put of the IF amplifier; the other 
ALSEP signals appear at 10 MHz 
+ ^0~^i)' ^ ~ ' 2, 3, 4. (See text.) 
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receiver, because its frequency conversion stages have 
"high" local oscillators, reverses the RF spectrum, so 
that, for example, the signal of an ALSEP whose S-band 
frequency is 1 MHz above the frequency of the reference 
ALSEP will appear at 9 MHz in the IF band. Table 3.1 shows 
the radio frequencies of the five ALSEPs, and the inter- 
mediate frequencies for two possible choices of the 
reference ALSEP. 


STDN 

ALSEP 

NUMBER 

APOLLO 

MISSION 

NUMBER 

S-BAND 

RADIO 

FREQUENCY 

USB RCVR 
(MHz) 

I.F. 



(MHz) 






A1 


A2 




REF 


REF 

A1 

12 

2278.5 

10.0 


9.5 

A2 

15 

2278.0 

10.5 


10.0 

A3 

16 

2276.0 

12.5 


12.0 

A4 

14 

2279.5 

9.0 


8.5 

A5 

17 

2275.5 

13.0 


12.5 


Table 3.1 Radio frequencies and USB receiver inter- 
mediate frequencies of ALSEPs for two 
possible choices of reference. These are 
nominal frequencies. Actual frequencies 
may deviate from the nominal by about 
20 kHz. 

Several other elements of the USB system are used by the Dif- 
ferential Doppler Receivers (DDRs) . The station primary frequency 
standard (PFS) , a cesium atomic-beam device, serves as a reference for 


two Hewlett-Packard model HP5100 frequency synthesizers, 
each of which provides a local-oscillator signal to one 
of the DDRs (see Section 3.2). The output of each DDR is 
in turn connected to one of the two input channels of the 
USB system Tracking Data Processor (TDP) . The heart of 
the TDP is a pair of digital counters which count, con- 
tinuously, cycles of the input signals. At selected time 
intervals controlled by the station clock, the readings 
of these counters are sampled non-destructively and recorded 
on paper tape. The TDP has two modes of operation. In 
the "single" mode only one channel is sampled and recorded, 
with the sampling interval selectable between 6 seconds and 
10 minutes. In the "dual" mode, the two channels are 
sampled alternately, and, in this case, the maximum 
sampling rate for each channel is once per 12 seconds. 

For ALSEP differential VLBI observations, the TDP is 
usually operated in the dual mode in order to count and 
record the output signal from two DDRs. 
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3.2 The M.I.T. Differential Doppler Receiver 

The Differential Doppler Receiver (DDR) was developed 
by Hinteregger and Counselman to enable the TDP to count 
directly the difference between the frequencies received 
from two ALSEP transmitters. Although it is possible to 
determine this difference without a DDR by subtracting 
frequency readings from two USB receivers with one USB 
receiver locked to each ALSEP signal, the results ob- 
tained in this way are much less accurate than those which 
can be obtained with the DDR ( compare Figures 1.3 and 1.4) . 

The DDR receives as input from the USB receiver the 
10 MHz reference signal, which is locked in phase with 
the reference ALSEP carrier signal, and the wideband 
(7-13 MHz) IF signal, which contains the carrier signals 
of all the ALSEPs whose frequencies are within approximate- 
ly ±3.3 MHz of the frequency of the reference ALSEP (see 
Section 3.1). In the DDR, one of these ALSEP signals 
(other than the reference) is selected and is mixed 
with a local-oscillator signal derived from the USB 
receiver's 10 MHz reference signal to obtain a beat- 
frequency signal which has the phase of the reference 
ALSEP carrier minus the phase of the selected ALSEP car- 


rier. The DDR's output, which is counted by the TDP 


comes from a voltage-controlled-multivibrator which is 
phase-locked to an integer multiple (360) of the beat- 
frequency signal. To explain more completely the opera- 
tion of the DDR, we refer to Figure 3.2. 


From HP5100 
Synthesizer 
f 


From USB 
Ref. Osc. 
10 MHz 


-r 16 


DIGITAL SELECTABLE- 
SIDEBAND MIXER 


Output to TDP, 


LOOP 

FILTER 


From USB 
7-13 MHz IF 


10 MHz + f s /16 


|LSB MIXER 

8 kHz 

fi-fo+V 16 ' 

LOW 

PASS 


PHASE 

DETECTOR 


Figure 3.2. Block diagram of the Differential Doppler 
Receiver (DDR) . In the 7-13 MHz wideband 
IF input from the USB receiver, the ALSEP 
carrier signal to be tracked differentially 
by the DDR appears at a frequency of 10 MHz 
- fi + fo- The HP5100 synthesizer frequency, 
fs» is set, and the sideband of the DDR's 
digital mixer is selected, so that the chosen 
ALSEP signal emerges from the DDR's LSB 
mixer between 2.5 and 5.5 kHz. The voltage- 
controlled -mult ivibrator (VCM) is phase- 
locked to this signal in a frequency-multi- 
plying (x 360) loop. 



Selection of the desired ALSEP signal for phase- 
locked tracking by a DDR is accomplished by a human 
operator who must first have measured, by means of a 
pair of USB receivers, the approximate difference be- 
tween the reference and the other ALSEP' s carrier fre- 
quencies. Based on this measurement, the operator sets 
the HP5100 synthesizer which is connected to the DDR, 
and he selects the appropriate sideband for the digital 
mixer which generates the DDR's local-oscillator (LO) 
signal, so that the desired ALSEP carrier signal will 
be converted within the DDR to approximately 4 kHz, 
near the center of the 2.5 to 5.5 kHz tracking range 
of the DDR's phase-locked loop. For example, if Al is 
the reference ALSEP and it is desired to lock a DDR 
to A2 , Table 3.1 shows that A2 ' s IF carrier frequency will 
be approximately 10.5 MHz, or 0.5 MHz on the upper side 
of the 10 MHz reference frequency. Referring to Figure 3.2 
we see that the frequency of the HP5100 synthesizer must 
be set in this case to 8.064 MHz and the upper (+) LO 
sideband selected, to obtain a DDR local-oscillator 
frequency of 10 MHz + (8.064 MHz/16) = 10 MHz + 0.504 MHz - 
10.504 MHz. The desired ALSEP carrier signal, at 10.5 MHz, 
is now 4 kHz lower than the LO frequency and will be 
converted to 4 kHz at the output of the final, lower- 
sideband (LSB) mixer. The output of this mixer is then 



low-pass filtered so that all noise signals farther than 
± 4 kHz from the chosen ALSEP’ s carrier frequency, and 
all other ALSEP signals, are rejected. 

If the ALSEP to be tracked by the DDR has an IF 
frequency on the lower side of 10 MHz (for example A4 
when the reference ALSEP is Al) , then the operator must 
switch the DDR’s digital SSB mixer to lower sideband (-) . 

The operator sets the synthesizer, in general, accord- 
ing to one of the following formulas: 

f = 16 (f n - f +4 kHz) Upper sideband case (3.2.1a) 

s 01 (f -l < V 

or 

f _ is /f - f - 4 kHz) Lower sideband case (3.2.1b) 

s 10 (f x > f 0 ) 

where f Q and f J are the measured S-band frequencies of 
reference and the other ALSEP , respectively. 

Setting the synthesizer and the sideband-selection 
switch are the only manual operations required in order 
to place the desired ALSEP carrier signal within the 
phase-lock tracking range of the DDR. However, it is 
usually also necessary for the operator to slew the VCM's 
frequency manually to the neighborhood of the proper 
locking point, 360 times the converted ALSEP carrier 
frequency, or nominally 360 x 4 kHz — 1440 kHz. 


A frequency counter is used by the operator to monitor the 
VCM frequency, which, once having been locked, will track 
any subsequent ALSEP frequency variations within the 
limits of 4 ± 1.5 kHz for the input to the phase detector, 

or 1440 ± 540 kHz for the VCM output. In general the DDR 

(VCM) output frequency, f D , in terms of the (constant) 
offset synthesizer frequency, f g , and the two (varying) 

ALSEP frequencies, is given by 

f = 360 (f. - f n + f /16 ) Upper sideband case (3.2.2a) 

f D = 360 (f 1 - f Q - f s /16 ) Lower sideband case (3.2.2b) 

We note that the frequencies of the local oscillators 
used in the USB receiver to convert the ALSEP signals from 
S-band (2.3 GHz) to the 10 ± 3.3 MHz band do not appear 
in this equation because the difference between ALSEP fre- 
quencies is not affected by these conversions as long 
as both ALSEP signals pass through the same mixers, etc. 

Also, the exact value of the "10 MHz" reference oscillator 
frequency, which is usually different from 10 MHz by 
about 50 to 100 Hz, does not matter: By virtue of the 

USB receiver's phase lock, the reference ALSEP' s IF 
carrier is locked to this reference oscillator, whatever 
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the latter’s frequency is; the first LO of the USB receiver is 
controlled to maintain this condition. If the frequency 
of the "10 MHz" reference oscillator is higher than 
10 MHz, so will the reference ALSEP's carrier be, and 
the other ALSEP's carrier frequency in the IF band will 
be increased by precisely the same amount . But in the 
DDR the second ALSEP's carrier frequency is subtracted 
from the frequency of the "10 MHz" reference oscillator, 
plus or minus f gyn /16. Thus, any variation in the "10 MHz" 
reference oscillator's frequency is cancelled out. The 
only frequencies remaining in Equation (3.2.2) are the 
ALSEP frequencies (their difference), and f syn / 16 « It 
is crucial, therefore, that the synthesizer be very 
stable, and be operating from the station's cesium-beam 
primary frequency standard. 

Since frequency is the time -derivative of phase, 
it follows that all of the above statements regarding 
insensitivity and sensitivity to the frequencies of Various 
oscillators apply with equal force to phase variations of 
the respective signals. The insensitivity of the DDR's 
output to phase variations within the USB receiver's 
LO system is the DDR's "secret of success". If we at- 
tempted to use separate USB receivers to measure the 


frequencies of two ALSEPs separately, in order to determine 
the frequency (or phase) difference, we would be susceptible 
to phase instabilities of both USB receivers. With a DDR, 
we are susceptible to no significant instrumental errors 
other than those of the DDR, the synthesizer, and prime 
frequency standard. In practice this has meant a reduction 
of two orders of magnitude in overall instrumental errors. 

The DDR is assembled from 15 integrated circuits, 

3 transistors, and associated passive circuit elements on 
a single board measuring 4 by 6 inches. The digital 
single sideband mixer is of novel design, using only 
digital logic devices to generate selectively either the 
sum or difference of the 10 MHz and offset frequency 
inputs [12] . The analog lower sideband converter is based 
on the design of A. E. E. Rogers [27] . The phase- 
locked loop is a conventional second-order type with 
lead-lag compensation, but using an exclusive-or logic 
gate as a phase detector. This type of phase detector, 
together with the use of differential circuitry in the 
loop filter, yields excellent performance with both 
strong and weak signals, with low sensitivity of the 
output signal phase to the input signal amplitude and to 
environmental factors. The divide— by— 16 in the input cir 


cuit for the offset frequency signal serves to reduce the 
effect of the phase drift of the synthesizer by that factor. 
With a 30-foot-diameter antenna and a 200°K system 
temperature (typical for an STDN station looking at the moon) 
the theoretically-calculated ALSEP carrier signal-to-neise 
ratio is -3 db in the 8-kHz bandwidth input to the phase 
detector, and +22 db in the +50 Hz closed-loop bandwidth of 
the phase-locked loop, which implies that the theoretical rms 
phase noise in the output of the loop should be 7°. The 
actually-measured loop performance is in good agreement with 
this calculated figure. Laboratory 

measurements that compared the outputs of two prototype 
DDRs have shown less than 1° of phase drift in 24 hours, 
including rapid temperature changes of 10°C to one 0OR. A 
field test in which two DDRs using separate synthesizers 
were locked to the same ALSEP pair showed approximately 1° 
of relative phase drift in six hours. 
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3 . 3 Observation Procedure 

The performance of the observations involves a co- 
operative effort by personnel at two or more STDN stations 
at the Network Operations Center (NOC) of Goddard Space 
Flight Center, and at MIT. The observations are directed 
from MIT* through NOC. Fifteen minutes before the begin- 
ning of each DDR "track", NOC calls the MIT Test Director 
and connects him to a voice communications network which 
includes each of the stations scheduled to participate on 
that day. At this time the Test Director tells the 
stations which ALSEP is to be used as reference and which 
ALSEP is to be tracked differentially on each of the two 
DDRs. At each station the operator then locks a USB re- 
ceiver to the reference ALSEP and connects the USB IF 
and 10 MHz reference signals to the DDRs. 

Before a DDR can be locked to its assigned ALSEP, 
the S-band frequencies must be measured for this ALSEP and 
the reference ALSEP, and the HP5100 frequency synthe- 
sizer setting calculated as described in Section 3.2. The 
two ALSEP frequency measurements are performed simultan- 
eously using two USB receivers with their associated fre- 
quency counters at the same station. As soon as the 

* 

From the initial installation of the DDRs in March 1973 
until July 1974, the period covered by the data included 
in this thesis, I directed most of the observations, as- 
sisted from time to time by Counselman and Hinteregger. 
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synthesizer frequency has been calculated, the (same) 
value is set into the synthesizers at both stations, and 
the operators lock the DDRs by slewing the VCMs to the 
anticipated locking frequency, usually 1440 kHz (see 
Section 3.2). Occasionally the DDRs will be found to lock 
at a frequency significantly different from that anti- 
cipated, due to drift of at least one ALSEP frequency 
between the time of the original measurement and the time 
the DDR is locked. In this case, a verification is 
always made that the lock is correct by again reading 
the two ALSEP S-band frequencies and the DDR output fre- 
quency and using Equation (3.2.2) to calculate the 
theoretical value of the DDR output frequency, which 
should agree with the measured value to within 3.6 kHz, 
i.e. to within 360 times the 10-Hz resolution of the USB 
frequency counters. If the ALSEP frequency drift is so 
rapid that an extrapolation indicates that the upper or 
the lower limit of the DDR's tracking range would be 
reached before the scheduled end of the observing period, 
the Test Director can instruct the station operators to 
change their frequency synthesizer settings and to re- 
lock the DDRs at a lower or higher frequency, in order to 
accommodate the drift. Alternatively, he can re-assign 
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that DDR to track a different ALSEP relative to the re- 
ference ALSEP. He may also choose a different ALSEP as 
reference. The frequency drift rate of an ALSEP is 
usually considerable only when the sunrise or sunset 
terminator is crossing the ALSEP, or for ~2 days there- 
after. Tracking assignments are therefore rotated to 
avoid tracking an ALSEP near its twice-monthly terminator- 
crossing times. 

After all of the stations have locked their DDRs and 
have begun recording data, the Test Director performs 
three additional checks to ensure that the data will be 
valid. First, he compares the stations' DDR output 
frequencies to see that the frequency for a given DDR 
at one station agrees with the frequency for the corresponding DDR 
at each of the other stations to within a few kHz. This 
comparison is particularly reassuring if the 
ALSEP transmitter frequency is drifting significantly, 
because the Test Director can then see that the DDR fre- 
quencies at all of the stations are changing in exactly 
the same way. As a second check, the Test Director per- 
forms a verification calculation similar to that which 
may have been done by one of the stations using Equation (3.2.2) , 
but in this case he obtains ALSEP S-band frequency readings 
from a station other than the one which performed the 
original calculation of synthesizer frequency. This check 
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is useful, for example, to determine that the DDRs are 
not locked on a sideband of the ALSEP carrier signal as 
a result of a mistake in the calculation of f s * Finally, 
an operator at each station reads to the Test Director 
two successive values of the TDP output on each channel. 

By subtracting successive N-counts and dividing by the 
time interval between them, the Test Director can deter- 
mine if the rate of increase of the recorded N-count matches 
the output frequency of the DDR which should be connected 
to that input channel of the TDP. The entire process of 
locking the DDRs and checking the output for two or three 
stations usually takes about one-half hour. Once this 
procedure has been completed, the MIT Test Director is 
disconnected from the communications network, but 
remains on call in case a problem arises later. The 
stations remain in contact throughout the observation 
period in order to compare DDR frequencies once each hour. 

3.4 Scheduling and Data Collection 

Since their inception in March of 1973, the ALSEP 
differential VL B I observations have been scheduled by 
Goddard Space Flight Center* on the basis of non-interfer- 
ence with other NASA tracking requirements. In- spite of 


Network support for the experiment wes^arran^ed God 
dard by I. M. Salzburg, J. W. Ryan, and H. W. Stone 
sifer, with the general support of Dr. J. M. ciarx. 
Director . 
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this restriction, observations have been obtained at 
least twice per week during most of this period. There 
have been several times when no observations were made for 
two- or three-week periods, but the only extended gap 
was from mid-April to the first of July in 1973 during 
the first Skylab mission. A complete list of observa- 
tions is given in Appendix A. 2. 

Because of the uncertainty of scheduling, our philo- 
sophy during the first year of the experiment was to ob- 
serve at any time that two or more stations could be 
scheduled for four or more hours. During this year, data 
were recorded on paper tape at the maximum TDP sampling 
rate, once per 12 seconds on each channel. The paper 
tapes were mailed by the stations directly to MIT where 
the data were transcribed onto magnetic tape for pro- 
cessing. 

By the spring of 1974, when a large volume of data 
had been collected and the frequency of the observation 
periods had become consistent (and none of the ALSEP 
transmitters had died) , we felt able to choose the ob- 
serving times and stations on a more scientific basis, fol 
lowing the geometrical analysis of Section 2.2. This 
analysis demonstrates the importance of performing ob- 
servations with both north-south and east-west baselines 
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to obtain maximum sensitivity to right ascension and 
declination separations, respectively, of the ALSEPs. 

The station pair of Madrid, Spain, and Ascension Island in 
the South Atlantic fits the north-south baseline criterion 
very well (see Figure 1.5). The best east-west station 
pair of those available would be Merritt Island, Florida, 
and Goldstone , California, if these stations could be 
scheduled for the entire period during which the moon 
is mutually visible, often ten hours or more. However, 
when other NASA tracking requirements limit the scheduled 
period to less than six hours, a baseline of Madrid or 
Ascension paired with either Merritt Island or Corpus 
Christi, Texas, is better by virtue of having a larger 
equatorial component [see Equation (2.2.7)]. In light 
of this analysis, and in an effort to optimize the expen- 
diture of resources by both the STDN and the MIT group, 
in May 1974 we reduced our schedule request to twd tracks 
per week, one using Ascension and Madrid, the other using 
the east-west baseline which was available for the long- 
est track during the week. The TDP sampling rate was 
also reduced to one observation per minute from each DDR. 
The lower sampling rate still provides enough samples 
during an observation period to ensure that short-period 
fluctuations will not be the dominant source of error 
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in the final parameter estimates. The reduction of the 
sampling rate also allowed the data to be transmitted 
from the stations via teletype to Goddard Space Flight 
Center and recorded there on magnetic tape, thus decreasing 
significantly the time and effort required for us to pre- 
the data for processing in PEP. 


pare 
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CHAPTER IV 

Data Processing Techniques 

4.1 Methods of Parameter Estimation 

The estimation of parameter values from the observa- 
tions is accomplished in PEP by weighted least-squares fitting. 

The theoretical basis for this method may be found in 
several texts (see [38], for example), and is not re- 
peated here. The reader is also referred to Ash [6 ] 
for a discussion of the implementation of this estima- 
tion scheme in PEP. In this section I will describe 
briefly the equations used and discuss their application 
to the data from this experiment. 

Suppose we have the vector of observations v , of 
dimension n, and a mathematical model which gives the 
theoretical value of the observations as a function of 
the parameter vector x, of dimension m(- n) . Then the 
vector of the observed values may be written as 
a sum of the theoretically computed vector y^ and the 
measurement-error vector e. i 

Yo = Y c (x) + £ (4.1.1) 

In weighted least-squares fitting, the estimate of the 
parameter vector, denoted by &, is the value of x that 
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minimizes the scalar quantity 

[v - v,(x)] T w [y_ - y_(x)] (4.1.2) 

=-o =-c — — — o — c — 

where W is the (n x n) weighting matrix. The diagonal 
_ 2 

elements of W are given by 1/6^, where 6^ is the assumed 
standard error in measurement i_, and the off-diagonal 
elements are zero (see below) . 

A closed-form solution for £ exists if the depen- 
dence of y on x is linear in x. If y is nonlinear, but 
=~c — — c 

well-enough behaved for small variations in x, the problem 
can be linearized. Then if a p riori values of the para- 
meters, denoted by x°, are available which are suffi- 
ciently close to the true values, we can write 

y = y ( x° ) + A ( x - x° ) + £ (4.1.3) 

=Lo =-c ' — ' — — — 

where A is the n x m matrix of partial derivatives with ele- 
ments a. . — I ay^/Sx - ] . If Equation (4.1.3) is valid 

ID 3 x° 

and the columns of A are linearly independent, then 
the least-squares estimate of x is given by 

x = x° + (A T W - Y^x 0 )] (4.1.4) 

If x° is not sufficiently close to £ at the outset, 
then one must iterate, replacing x° with the right 
side of Equation (4.1.4), then re-evaluating that 


side. The solution is said to be converged when sub- 
stitution of the estimated values of the parameters back 
into the model produces no significant further decrease 
in the sum of weighted squared residuals defined by 
Equation (4.1.2). It is possible, in principle, that the 
solution so obtained will correspond to a relative, rather 
than to the absolute, minimum of Equation (4.1.2). In 
practice the parameter values are sufficiently well known 
a Priori that the solution will converge to the proper 
values . 

The statistical properties of the weighted-least- 
squares estimate depend mainly on the statistics of the 
errors of observation, e_. In particular, if the 
elements of e are samples of zero-mean, Gaussian, inde- 
pendent random variables, and if each observation is 
weighted inversely with its associated error variance, 
then the weighted-least-squares estimate is the maximum- 
likelihood estimate. In this case the covariance matrix 
of the error in £ is given by 

P= (A T W A)" 1 (4.1.5) 

If, further, we denote the elements of P by P^j/ the 
so-called formal standard error of the estimate of 
parameter i is 
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a i = *^ii (4.1.6) 

and the correlation coefficient between the errors in 
the estimates of and is 




(4.1.7) 


In practice, cs^ in Equation (4.1.6) can usually be considered 
only a lower bound on the uncertainties in the para- 
meter estimates. The actual observation errors are rrost often 
dominated by processes with significant systematic com- 
ponents. The effect of such systematic errors on the 
parameter estimates is discussed in Section 4.5. 

One source of significant measurement-error correla- 
tions in most counted-cycle observations has been elim- 
inated in this experiment by the way in which the ob- 
served quantity has been defined. It has been con- 
ventional in Doppler tracking to define the observable 
as the difference between successive counter readings 
divided by the time interval between them. The result- 
ant from this operation is the average frequency over the 
interval, a quantity which can be interpreted in coherent 
Doppler tracking as the average (two-way) Doppler shift 
of the signal between transmission and reception. How- 
ever, because the error in a given sample of the counter 


contributes in an equal and opposite manner to each of 
two neighboring average Doppler values, the errors in 
successive measurements will have a correlation of -0.5 
if other noise sources are ignored. Using the simple three- 
parameter model described in Section 2 . 2 of this thesis, 
Shapiro et al. [32] have shown that the neglect of this 
correlation leads to a significant overestimate of the 
error in lunar baseline determination. The formal 
standard errors obtained from a single six-hour obser- 
vation period will be too large by a factor of n, where 
n is the number of observations, and the a ctua l, error 
standard deviation will be too large by a factor of /n. 

If the observable is chosen to be the accumulated N- 
count, formed by subtracting the counter reading at the 
beginning of the observation period from each successive 
reading, this problem of correlated sampling errors can 
be easily avoided (see Section 4.2.1). The -0.5 cor- 
relation between the error in each measurement and the 
error in the initial counter reading is eliminated by 
solving for an unknown constant for each series of unin- 
terrupted observations. This "bias" parameter will ab- 
sorb a number of errors which are constant throughout a 
tracking period (see Section 4 . 5 ). 
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4 . 2 Theoretical Model for the Observable 
4.2.1 Basic Formulation 

By counting cycles of the output signal of 
the voltage-controlled multivibrator (VCM) of the DDR, 
the tracking data processor (TDP) at each station measures 
the phase of this signal as a function of time. As describ- 
ed in Section 3.2, the change in the phase over a given 
time interval is just 360 times the change in the phase 
difference between the carrier signals received from the 
two ALSEPs, plus a known bias which is introduced by the 
frequency synthesizer. Some random measurement noise is 
also present. Since the initial reading of the TDP is 
arbitrary, without loss of generality we may define this 
reading to be zero. Then the reading of the TDP at 
station j at time tj is 


- +rlj (t j ) + Vj'V + Vj (t j )] 

(4.2.1) 


where <P . .(t.) is the phase of the signal received 

from the i th ALSEP at the 3 th station at time t ^ ; 

A . is the known bias; and <j> M . is the random noise. Note 
Y Bj N J 

that, by definition Nj(t°) = 0, where tj is the time of 


i 

sr • 


f 
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the initial reading of the TDP. The TDP counts whole 
cycles of the signal from the VCM, so that in Equa- 
tion (4.2.1) the unit of Nj is l/360th of a cycle, or 
1° of phase. The error due to the integer-degree trun- 
cation is included in the noise term, 4^ j • 

The times in Equation (4.2.1) are the coordinate 
times, rather than the station clock readings, when the 
cycle-count is sampled by the TDP. The relationship 
between coordinate time and the station time is given 
in Appendix^ .and discussed in some detail by Ash [ 6] . 
Briefly, coordinate time is the independent variable in 
the equations of motion which we integrate to obtain the 
positions and velocities of the moon, the earth, and 
the planets. Coordinate time and atomic time differ by 
small diurnal, monthly, and annual periodic terms. The 
station clocks are set, within 10 to 20 ysec, to Co- 
ordinated Universal Time (UTC) , which has a well-known (defined) 

relationship with atomic time. 

The signal received at the time t ^ was transmitted 

by the i th ALSEP at the time t ± given by 

(4.2.3) 

fc i “ ' T ij 

where t. . is the propagation time delay. The phase 
of the received signal, 4> ri j/ 1# Related to the 
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phase of the transmitted Signal, 4> ti » by 


♦rij'V “ ' 

Thus, the TDP reading at 
in terms of the transmitted phases : 

W = t* t 2 (t j ‘ T 2j ) " ♦tl^j * T lj ) 

+ 4»_. . (t . ) + 4> n . (t . ) 1 - [similar terms 
Dj J J J 

for t?1 

The differenced N-count observable * is defined 
as the simple arithmetic difference between the TDP 
readings, or N-counts , at two stations, taken at 
identical readings of the two station clocks. Suppose 


x .) (4.2.4) 

ID 

station j may be expressed 


This observable might equally well have been terme 
the "doubly-differenced N-count", since the phase is 
differenced both between stations and between ALSEPs. 

I have chosen the simpler term inasmuch as the difference 
between ALSEPs is taken in the DDR before the counting 
is performed. 


-69- 


that when these clock readings both equal T™ , the coordinate 

times are t™ and t^ at station j and station k, res- 

3 .. . th 

pectively. Then the differenced N-count at this m 

sample time" is 

AN^ k = Njttj) - N k (t k ) (4.2.6) 

Substituting Equation (4.2.5) into Equation (4.2.6), we obtain 

an^ = u, t 2 (t- - t^) - * t 2 (t; - x; k )i 
"[♦ tl (t? ' ^j) - *tl<< - 
+ [(f, Nj " ^Nk 31 “ {similar termS 


for m = 0} 


(4.2.7) 


where the superscript m on and 4> N j denotes that the 

coordinate time argument is t™. The 4> B j and ^Bk " bias 
terms are omitted deliberately. These terms cancel since 
the synthesizers that determine the bias frequencies at 
both stations have the same setting : the synthesizer 

and the clock at a given station use the same reference 
frequency standard. 


We assume that the transmitter frequency as a 
function of time can be represented by 


4> ti (t + At) = 4> ti (t) + i ti (t)At 


+ | $ ti (t) At 2 + . . . 


(4.2.8) 


where At is some as yet unspecified interval. Using 
Equation (4.2.8) we expand Equation (4.2.7): 

•.jn r p i / . iti in v i / 1 m m \ t < . in m * / 1 m , m m . m » 

*jk = ' ' , ’t 2 (t j -T 2 j' - ' t 't 2 <t j- T 2 j )( V t j- t 2 k +T 2 j ) 


1 j / 1 m m » / 1 Tn , m m . m \ 2 i 

- J 4’t2 (t j' T 23 ) (t k _t j ~ T 2k + T 2j “ 


,.m m 


, m _m 


m m . m jrci m , m 


r i / , iu xu . , / . xu xu » I / . xii xii * /j_*“ . ui . in \ 

- c<f> tl (tj - T ij>- +ti (t j- T ij ) - ,f 'ti (t j- T ij )< V t j- T ik +T ij > 


or 


1 1 / 1 m m y /x.^ i_m ,m .2 i 

" 2% (t j “ T 1 j > (t k _t j" T lk +T lj > “ * * * 1 

+ ” {similar terms for m = 0} 


(4.2.9) 


AN jk 


ri / , m m . / m m , i < . m _«» » < _m \ 

= Ty) (ty-T^) - * t2 (t j IjI 


m m . , m m 


jn m 


m m \ i / « m . nil 


* r I # *xil \ 1 xii y i /jin 

+ [4.ti(tj-Tij)-* t 2 <t j-T 2j ) ] (Vtj) 
lv m m v m m m . m .2 

- 2 ^ t 2 (t j- T 2 j ) <V t j- T 2 k +T 2 j> 

. '1 jl ,.m m . ,,m .m m . m ,2 , 

* 2 1 k (t j > <' W T ik +T ij > + • • • 


+ [£>-&] ) ” {similar terms for m = 0} 


(4.2.10) 
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The terms in Equation (4.2.10) have been listed in 
descending order of magnitude. The geometric information 
in the observable is contained primarily in the terms m 
the first line. The main computational task in the 
determination of the theoretical value of the observable 
is the evaluation of the propagation delays in the ex- 
pressions of line 1. Subsection 4.2.2 deals with this 
evaluation. A second task, much simpler but equally im- 
portant, is to determine the transmitter frequencies, 
^ti ( t m -T^j) and * t 2 (t^j), which are used to convert the calculated 
delays to phase. Subsection 4.2.3 deals with the evaluation of 
these quantities. 

The terms in line 2 of Equation (4.2.10) account 
for the non-simultaneity of the receive times. The only 
effect modeled here is the approximately 1 ysec diurnal 
variation of the atomic-to-coordinate time difference. 

We have not attempted to account for cgnstant errors in 
clock synchronization, or for offsets of the rates of the 
clocks at the stations from UTC. The errors introduced 
into our model of the observable by the neglect of these 
effects are discussed in Section 4.5. 

The third and fourth lines of Equation (4.2.10) 
contain terms involving variations in the transmitter 


. . ....... 






f requenc is s • Ws shall see m ths discussion of transmitter 
frequency characteristics in Subsection 4.2.3 that the 
effect of these terras is insignificant. 

Finally, the random measurement noise term is not 
included in the theoretical model of the observable 
except for the initial observation (m = 0) . For a given 
series of observations, the same unknown (solved-for) 
constant is added to the value of ANjj^ for all m — 1, 
to represent the constant contribution to AN^ of the 
error in the {m = 0} term in Equation (4.2.10), as dis- 
cussed in Section 4.1. 


4.2.2 Evaluation of Propagation Delays 

Each of the propagation delays, x™j , in Equa- 
tion (4.2.10) can be considered, with sufficient ac- 
curacy, to be the sum of a geometric (free-space) delay 
plus a small additional delay representing the effect 
of the propagation medium on the signal. The geometric 
delay is the magnitude of the vector 


1 

c 






+ r ti ( ti ) 


r e (t 3> ' 



(t. 

3 


) 


(4.2.11) 


where 


c is the speed of light; 

r m (t^) is the position of the moon's center of mass 
with respect to the solar-system barycenter 
at the time of transmission t^; 
r ^(t^) is the position of transmitter i_ with 

respect to the moon's center of mass* at t^; 
r e (tj) is the position of the earth's center of mass 
with respect to the solar- system barycenter 
at the time of reception t ^ ; 

and 



is the position of the observing site with 
respect to the earth's center of mass* at 



The solar - system barycenter is used as the origin of the 
coordinate system used for the calculation of R^j so that 
our coordinate system may be indistinguishable from an 
inertial system. The directions of the coordinate axes 
of the system are defined in theory by the mean equator 


Strictly speaking the vectors r . . and r . 

l1 5 J 


should be ex- 


pressed properly in the solar-system-barycentric system, 
I ignored this refinement and added the coordinates 
(even though expressed in the body-fixed systems) 
directly in Equation (4.2.11). The resulting errors 
should be completely negligible. 


-74- 


and equinox of 1950.0. Operationally, the definition 
of the coordinate system is provided by the models used 
for the orbits of the planets and the moon and the motions 
of the earth and moon about their respective centers of 

maSS * ThS planetar y ephemerides which I used were gen- 
erated at MIT by adjusting the system of planetary masses 

and the initial conditions of a numerical integration of the 
Planetary orbits to fit ^y years of optical , rato/ ^ 

tracking data [2, 4, 5] . The particular ephemeras tape em- 
ployed in processing the data for this thesis is named 
NBODY311 , dated 16 August 1969. My lunar ephemeris re- 
from fitting a numerical integration to four and 
one-half years of laser ranging data [ 1 , 33 , 34 ] , using 
the equations originally developed at MIT by M. a. slade 

The model for the motion of a transmitter on the 

lunar surface with respect to the moon's center of mass 

included both the mean rotation and the physical libra- 

tion. This model is of particular interest for this thesis 

and is described in detail in Section 4.3. Finally, the 

calculation of the positions of the observing sites in 

the earth-centered 1950.0 coordinate system, including 

precession, nutation, variation of UT 1 an H n i 

1 U1 * 1 » and polar motion, 

has been documented by Ash [ 6 ] . 
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Given the coordinate time, t, of reception of a signal. 


the geometric delay (t) is determined iteratively: 


. _ lim x *}^) (t) 

T ij ' k->°° lj 


( 4 . 2 . 12 ) 


where 


T ij k+1) (t) 




t! i k) (t), t) 

1 j 


k = 0 , 1, 


and 


r! (t) = 1.25 seconds 
ID 


(4.2.13) 

(4.2.14) 


The additional propagation delay due to the at- 
mosphere was not included in the light-time iteration 
but was added as a small correction to (t) to yield 
t (t) The amount of propagation delay due to the neutral atmosphere 

ij • 

was calculated for each site from the monthly-average data 
for the zenith delay compiled at the Goddard Space Flight Center 
[28 ] , using a model developed by Chao [10 ] for the effect 
of elevation angle. No attempt was made to model the 
effect of the ionosphere on the propagation delays. Es- 
timates of the errors from both the neutral atmosphere 
and the ionosphere are given in Section 4.5. 

4.2.3 Determination of Transmitter Freque ncies 

For the evaluation of the theoretical value 
of the observable by means of Equation (4.2.10) we also 
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re quire the values of the transmitter frequencies and their 
rates of change. The accuracy with which these frequen- 
cies must be known, for a given uncertainty xn the cal- 
culated value of AN^j , is set by the magnitude of the 
propagation-delay difference (t^ - T ±k ) in the first line 
of Equation (4.2.10). This difference depends on the 
rotation of the earth, the positions of the stations, 
etc., but is usually on the order of a few msec and is 
rarely greater than 10 msec. For the theoretical ob- 
servable to have an uncertainty of less than 1 , or 0.003 
cycles, the transmitted frequency must be known to within 
(0.003 cycles) * (.01 sec) = 0.3 Hz, or to about 1 part 

in 10 10 . 

The ALSEP frequencies are not known a priori with 
nearly this accuracy. These frequencies can vary un- 
predictably by several kHz from one day to the next, 
by several hundred Hz to about 1 kHz over 1 hour, and 
by about 1 Hz in 1 sec, representing even in the last 
instance a change of 5 parts in 10 10 . Thus it is neces- 
sary to monitor the frequencies of the signals received 
from the ALSEPs at one (or more) of the tracking 
stations, and to calculate the transmitted frequencies 
from the measurements of the received frequencies by 
accounting for the Doppler shift between the transmitter 


and the receiver. How accurate can this procedure be? 

The accuracy of the cesium-beam reference frequency 
standard at the station exceeds our requirement of 1 
part in 10 10 . The uncertainties of the a priori posi- 
tions of the tracking stations, the geocentric ephemeris 
of the moon, etc., which are needed to calculate the 
Doppler shift, are smaller than the equivalent of 1 
second of arc, which corresponds to an uncertainty in 
the determination of the transmitter frequency of less 
than 1 part in 10 11 . Evidently, then, there are no 
^.¥P^ amen tal difficulties involved in using measurements 
of the frequency of the signal received at a tracking 
station as a function of time to determine the transmitted 
frequency as a function of time , with the accuracy desired 
for the calculation of Equation (4.2.10). However, a 
significant practical difficulty for the observations 
analyzed in this thesis has been that the only available 
received-frequency measurements have been relatively crude. 

At one tracking station, the operator uses a digital 
counter to measure the received frequency of the refer- 
ence ALSEP approximately every 15 minutes. The uncer- 
tainty in these measurements, including the uncertainty 
associated with interpolation of the frequency within the 
15-minute intervals between measurements, is about 100 Hz. 
From the foregoing discussion it would 


appear that the 


-78- 


uncertainty of 100 Hz in this frequency measurement would 
lead to an uncertainty of a full cycle, or 360° of phase, 
in the evaluation of Equation (4.2.10). Fortunately, 
a much more accurate measurement is available of the 
difference between the received frequencies of the two 
ALSEPs — from the output of the DDR. As long as the 
difference is known accurately, a relatively large common 
error in the two frequencies can be harmless, because the 
error introduced in the term in Equation (4.2.10) in- 
volving <j> t2 is very nearly cancelled by the error intro- 
duced in the term involving 4> tl * In fact ' since 
generally 

|(t»j - t^) - (T™j - x^ k ) i * 50 ysec (4.2. 

a simple common-mode frequency error of 100 Hz results 
in less than a 2° error in the evaluation of 
Equation (4.2.10). But a much smaller error in the 
measurement of the frequency difference can be much 
more serious because no comparable cancellation occurs. 
How well is the difference between the received fre- 
quencies known? For the observations analyzed in this 
thesis, the uncertainty associated with the difference- 
frequency measurement was of the order of 1 to 10 Hz, 
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due mainly to the fact that it was not possible to record 
measurements of the difference frequency with good time- 
resolution. The phase of the output of the DDR, which 
alone contains the needed information, was usually sampled 
and recorded by the TDP only at 60-second intervals. In 
fact, for most of the observations only every third sample 
was actually processed, in order to conserve computer 
time. Thus the frequency-measurement interval in most 
cases was 180 sec long. For historical reasons it was most 
convenient to use the difference-received-frequency 
measurement derived from the interval immediately preceding 
t™, to remove the a priori Doppler shift from the received 
frequency, and to use the resulting value of transmitted 
frequency, assumed constant , in the evaluation of Equa- 
tion (4.2.10). Usually the ALSEP transmitters were suf- 
ficiently stable that, even though the frequency-measurement 
time interval was centered ~90 sec before tj 1 , the error 
in Equation (4.2.10) due to the use of an inexact value of 
transmitter frequency was less than 5° (rms) . However, in 
some cases the drift of the transmitter frequency was 
rapid and systematic, with a drift rate of ~0.5 Hz sec ^ 
persisting over several minutes or more . In such a case 
the just-described procedure would yield a frequency which 
was incorrect by ~45 Hz, and an error in Equation (4.2.10) 
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possibly as great as (~45 Hz) • (360°/cycle) • (0.01 sec) 
a 160° , which would be intolerable. In cases of such 
rapid drift, I usually processed every data point in- 
stead of every third one so that the frequency-measure- 
ment time interval was centered only 30 sec before t™; I 
also used a slightly more complicated algorithm to deter- 
mine the transmitted frequencies, based on the use of 
three successive TDP readings to estimate the drift rate 


of the transmitted frequency. 

In all cases, I programmed PEP to print out the rates 
of change of the transmitted frequencies and the calculated 
effect of those rates on the observable. If the error in a series of 
observations due to erratic behavior of the transmitter seemed likely to 
exceed a few degrees of phase, the observations were 
deleted from the later fits. 

The specific algorithms used for the calculation of 
transmitter frequencies are as follows. We calculate the 
average transmitted frequency of the reference ALSEP over 
the transmit-time interval corresponding to the receive 
time interval t™ 1 to t™, by 


, (n 4 ) 

tl 


(m-t) 

f i • 

rlD 


m— ; 


<1 - By ) 


(4.2.16) 


where f rlj 


is the received frequency of the reference 


ALSEP at the midpoint of the receive-time interval, deter- 


mined by linear interpolation from the digital counter 
readings recorded by the operator at 15-minute intervals, 

and 

_m-l 

(4.2.17) 


1 m 

ft m ~2 = T ij ~ ij 
^ij , m ,m-l 


, ill .ill 

" D 

is the average first-order Doppler shift over the inter- 
val t m-1 to t m . A more precise calculation of 


(m-i) , _ . . 

f 2 is not warranted since the accuracy of the 

result is limited at present by the accuracy with 

which the received frequency can be determined for the 

reference ALSEP . 

As indicated previously, the transmitted frequency 
of the second ALSEP is calculated from the frequency of 
the reference ALSEP using the DDR output, so that the 

large error in f fcl due to the crude measurement of 


< m 4> 


(m-i) 


f 2 is present equally in f A • From Equation (3.2.2), 
rlj ^ 

the average frequency of the DDR output at station j_ 
over the interval t m_1 to t m is given by 


T--' 
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( m-|) A yt^) - y fj’ 1 ) 


Dj 


m _ .m-1 


= 360 [f 


(m-|) 

r2j 


f 


(m-|) 

rlj 


+ f . 1 

b: j 


(4.2.18) 


where f . is the bias frequency which can be either 
B j 

positive or negative depending on the setting of the 
digital sideband selection switch in the DDR (see 
Section 3.2). Rearranging terms, we obtain an expression 
for the average received frequency of the second ALSEP : 


1 x 

T {m ~2 ] _ - f (m “2 ) f Dj _ f 
f r2j ” t rlj 360 Bj 


(4.2.19) 


transmitted frequency is 


(4.2.20) 


where 3™. 2 is defined by Equation (4.2.17). 

^ D 

The rates of change of the transmitter frequencies 

_(m-~) ^ _(m-|) 

are determined from the changes in f fcl ana f fc2 


From this value, the average 
calculated using 


"t2 


_(m-i) 

F r2j 

( 1 - 3 ^) 


between successive TDP sampling intervals. It is assumed 
that the rate of change of frequency is constant over 


that 
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two adjacent sampling intervals. Thus the transmitter 
frequency drift rate is calculated from 

/ i \ (m-i) _Jm-f) 

= (f - f I (4.2.21) 

ti l l(t m. t m-l ) + (t m-x_ t m-2 n 

If this rate is significant/ as discussed previously, 
then the transmitted frequencies used for the cal- 
culation of the theoretical observable in Equation (4.2.10) 
are extrapolated to correspond to the tin© of reception, 
t m , rather than to the midpoint of the previous sampling 
interval. That is, we use 

m + lj "m _ t n>-l) (4.2.22) 

r ti ti 2 ti' 

In any case, the value of f^-v, expressed in degrees 
of phase per unit time, is substituted for (tj - T ij) 
in the evaluation of Equation (4.2.10). The terms involving 
$ . and higher derivatives in Equation (4.2.10) are ig- 
nored . 


/ ' 





-84- 


4.3 Model for the Rotation of the Moon 

The dynamical equations of motion for the moon's ro- 
tation about its center of mass are the well-known Euler 
equations for rigid-body rotation: 

Au)^ + (C-B) ^ 2^2 = T 1 

Bo )2 “ (C— A) = (4.3.1) 

t (B— A) = ^3 


in which u^, w 3 , are the angular velocities of the 

moon about its principal axes; A^ B£ C are the moments 
of inertia about these axes; and T^ , T 2 ? and are 

- 4 - 

the corresponding components of the external torque, T, 
exerted by the gravitational fields of the earth and 
sun on the asymmetries of the moon's figure. The rotation 
i. usually parameterized by writing Equations (4.3.1) in 
the form 


+ X T 1 

“2 = 3w l w 3 + k T 2 

w 3 = -yu)-^ + qT 3 


(4.3.2) 


in which a = (C-B)/A, 3 - (C-A)/B, and y - (B-A)/C. 
These three ratios are not independent, but 
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a _ _ y)/(l - gy) . The torques are functions of 

g and y, and also of coefficients of the moon's gravi- 
tational potential above the second degree (see below) . 

In principle one could generate a model of the lunar 
rotation by directly integrating Equations (4.3.2), ex- 
pressing, with three Euler angles, the orientation of the 
moon's principal axes with respect to a set of axes fixed 
in inertial space (for example, the axes of the 1950.0 
equator-equinox system) . In fact, from the standpoint of 
modern numerical methods, this straightforward approach 
has many advantages. All models which have been developed 
thus far, however, introduce intermediate quantities in or 
der to take advantage of the fact that the moon's rotation 
mav be characterized by small-amplitude libra tions about 
a quasi-uniform "mean” rotation. This mean rotation was 
first described in 1693 by Cassini [9 ] , in terms of two 
"laws" : 

(i) The moon rotates uniformly about its axis of 
maximum moment of inertia with a rotational angular velocity 
equal to its mean orbital angular velocity around the earth. 

(ii) The lunar equator maintains a constant incli- 
nation with respect to the ecliptic; the planes of tne 
lunar equator, lunar orbit and the ecliptic intersect 
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along a common line, and the plane of the ecliptic lies 
always in between the planes of the equator and the orbit 
As a consequence of these "laws", it has been cus- 
tomary to describe the moon's rotation with respect to 
the rotating set of axes defined by the so-called mean 


ecliptic and mean equinox of date: the conventionally 
defined mean-of-date (MOD) ecliptic coordinate system (see 
below), usually the orientations of the moon's principal 
axes of inertia with respect to the MOD ecliptic coordinate 
axes are specified by three Euler angles, e, and as 
shown in Figure 4.1. The angular velocities v ^ , and 

u> 3 are related to the time-derivatives of these Euler 
angles by 

“i = sinesin <p - dcostp 

w 2 = si n0cos<J) + 0sin<t> 3 3 j 


w 3 - ^ cos0 + (J) 
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where Q is the longitude of the ascending node of the 

moon's mean orbit on the mean ecliptic of date, measured 

from the mean equinox of date along the ecliptic; I is 

the inclination of the lunar equator with respect to 

the irean ecliptic of date; ([ is the mean orbital longitude of 

the moon measured from the mean equinox of date; tt = 3.14159. . . ; a, 


p. and x are the angles describing the physical librations 
in node, latitude, and longitude, respectively. The angles 
fi and (£ are obtained conventionally from the Improved Lunar 
Ephemeris (see [17]). 


fi 

a 


259° 10 1 59779 - 5 r 134 °08 ' 31’.' 23T + 7'.'48T 2 
+ 0 008T 3 

270°26 , 02'.' 99 + 1336 r 307 °52 ' 59V 31T 
- 4'.'08T 2 + 0V0068T 3 


(4.3.5a) 

(4.3.5b) 


where T is measured in Julian centuries of 36525 

ephemeris days from Julian Day 2415020.0 (1900 Jan 0.5 E.T.) 

The libration model used for this thesis resulted 
from a numerical integration, performed by VI. S. Sin- 
clair, M. A. Slade, and J. G. Williams of the Jet Pro- 
pulsion Laboratory, of the equations of motion for the 
quantities la, p, and t. These equations were derived 
from Equation (4.3.2) by means of Equations (4.3.3) and 
(4.3.4), with the torque given by 
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T = M e r e x (VU) e + M s ? g x {VU) g (4.3.6) 

where M , M are the masses of the earth and sun, and 
0 S 

r g are the positions of the earth and sun relative to the 
moon, and (VU) is the gradient of the moon's gravitational 
potential at the center of mass of body x. The potential is represented 
by an expansion in spherical harmonics : 


U 


k 2 M 

a - l 

r n=2 


R m n 

j (— ) 

n 'r 


p (sinL) 

n 


+ I I [C nh c°s h* + S nh sin h«.]P nh (sinL)} 

n=2 h=l 

(4.3.7) 


where r, L, and £ are, respectively, the radius, latitude, 

V 2 

and longitude referred to the principal-axis system; 
is the gravitational constant (see Appendix A.l); M m is 
the mass of the moon; R m is an adopted value for the 
'moon’s equatorial radius; J n are the zonal-harmonic co- 
efficients; C nh and S nh are the tesseral cosine and sine 
harmonic coefficients, respectively; and P n and P nh are 
the Legendre polynomials and generalized Legendre 

functions, respectively. 

The gravitational harmonic coefficients enter 
Equation (4.3.2) multiplied by factors of and 

divided by A, B, or C. If the value of either J 2 or C 22 
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is known independently, e.g. from observations of spacecraft orbiting the 
noon, then this value may be combined with the values of 8 and y to 
determine C/M m R^ . Thus, the measurements of the third- 
and higher-degree harmonic coefficients derived from the 
physical libration will be uncertain by a common scale 
factor, with the uncertainty in scale set by the uncer- 
tainty in the independent determination of J 2 or ^22 or ' 
equivalently, in C/M^- The relationships among these 
quantities and their uncertainties are discussed by Williams 
et al. 140]. In the JPL libration model, J 2 has been 
fixed at a value of 2.03822 x lO -4 , determined from an 
analysis of Lunar Orbiter data. The adjustable parameters 

in the model are 6, Y / J 3 / C 31# C 32 ' C 33' S 31' S 32' and 
S 33 . Values for C 22 and C/jy^ may be derived from the 

adopted value of J 2 and from the values for 8 and y 
determined from the observations. A description of the 
procedure used at JPL to generate the integrated libration 
model is also given in [ 40] . The numerical values of 
the parameters in the particular model I have used are 
listed in Appendix A.l. 

The values of t, p, and o determined from the libra 
tion model, in combination with the conventionally defined 
values of <[ and define the relationship between the 
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selenodetic coordinates of an ALSEP , X, 
responding MOD ecliptic coordinates, X' 


and the cor- 
[see Figure 4.1 


and Equations (4.3.4)]. 


X = U X' 

where 


(4.3.8) 



cos^cos (f)-sinijjcos0sin<t> sin^cos(})4cosiJjcos0sin<J> -sinOsincf) 
-cosijjsin<|)-s.iji^cosecos4) -sini{Jsin^+cos4>cosecos4) -sin0cos4> 
-sint|jsin8 cosipsin© cx>s6 


(4.3.9) 


The transformation from MOD ecliptic coordinates 
to 1950.0 coordinates, x, is defined by Ash [3 ]: 


x = A T X* 


(4.3.10) 


where 


A = 


11 

p 21 COSe o +P 31 sine o 

P 31 cose o -p 21 sine | 

21 

p 22 COSe o +P 32 Sine o 

p 32 COSe o" P 22 sine ' 

31 

p 23 COSG o +p 33 sinE o 

P 33 cose o “P23sine 


(4.3. 11) 


In Equation (4.3.11), e Q is the conventional value for the 
mean obliquity of the earth's equator to the ecliptic, 
given by [17 ] 


e = 23°27’ 08726 - 467845T - 070059T' 
+ 07 00 18 IT 3 


(4.3.12 ) 


and are the elements of the precession matrix. The 

precession matrix used in PEP is derived from the con- 
ventional representation [17]; in addition, the matrix 
includes the possible effects of small additional angu- 
lar rates about the three axes of the 1950.0 system. These 
rates may be adjusted in the least-squares solution. For 
a complete definition of the precession matrix, see Ash 


[ 6 1 • 


Using the fact that U is orthogonal (hence, U 


= U T ) , we can combine Equations (4.3.8) and (4.310) to 
give the complete transformation from the selenodetic 
system to the 1950.0 mean equator-equinox system: 


x(t) = A T (t)U T (t)X 


(4.3.13 ) 


where the time dependence of x. A, and U has been 
denoted explicitly. 
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4.4 Partial Derivatives of the Observable with 
Respect to Parameters 

Using Equations (4.2.10) and (4.2.22), we can write 


the differenced N-count observable as 


^ - r t2 (T 2j - X 2k> 
- [similar terms for m = 0] 


(4.4.1) 


Changes in the values of the parameters in our model 
of the earth-moon system affect the propagation^delays , 
t • • (t m ) . If we ignore the weak dependence of f^ and 

f 1 " on the delays, the partial derivative of AN^ k 
with respect to a parameter, p, may be written 


3 AN 1 ?. m 

J* - f 
3p tl 


3x 


m 


9P 


3x 


m 

lk 


3p 


] 



[!^ti 

1 3p 


3x 


m 

2k 


3p 



- [similar terms for m = 0] 


(4.4.2) 


From the discussion above Equation (4.2.11) , 


3 T ij _ i 3 

3p c 3P 


(4,4.3) 


where |R ± j (t i , t j ) | is the magnitude of the vector 
pointing from receiving site ± to transmitter i on the 
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moon. We will denote the Cartesian coordinates of 
^ij(ti'tj) in the 1950.0 system by the mathematical 


vector {3x1 column matrix of the coordinates) x. .. 

-ID 

The unit of x^j is the light-second, the fundamental dis- 
tance unit of the observations. Then Equation (4.4.3) 
can be written 


3t. . 

hi 

3p 


9 x . 

— l 


D 1 _ 


3p 


T 
x. . 
- r 1 : 


SiJ 


T 
x. . 
-ID 

'SiJ 


9x. . 
-ID 

9p 


9x 

+ 

1 8P 


9x 


ti 


9p 


9x 

-e 

9p 


B ~sj 

9p 


(4.4.4) 


where x^, x^, x c t and x g j give the 1950.0 Cartesian 
coordinates of the vectors r jn (t i ) / ? ti (t i ), r Q (tj ) , and 
r . (t.) , which represent the positions of the moon, 

ALSEP transmitter, earth, and observing site, respectively, 
as defined in Section 4.2.2. 

4,4.1 Partials with Respect to Parameters of the 
Lunar Orbit 

ALSEP differential VLBI observations are 
relatively insensitive to the orbital motion of the moon. 
However, partials of the theoretical differenced N- 
count observable with respect to the 6 initial con- 
ditions of the lunar orbit and the mass of the earth- 
moon system were calculated , in addition to the par- 
tials of the theoretical laser ranging observable. The 
latter observable is, of course. 


t -i 
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highly sensitive to the orbital motion. For any one of the 

8 -ij 

6 orbital and 1 mass parameter (s) , the derivative g p 
in Equation (4.4.4) is simply 


9X . . rj 

-TAl = |_(x - x ) 

3p 9p -m -e 


(4.4.5) 


I integrated these partials 

numerically, simultaneously with the orbit of the moon, 
using the equations described by Ash [ 3 ] and Slade [36 ] 
The components of the partials at a particular epoch 
are determined by interpolation from values written on 

the integration tape. 

4.4.2 Partials with Respe ct to the ALSEP 
~ Selenodetic Coordinated 

From Equation (4.3.13) the selenocentric 

coordinates of an ALSEP transmitter expressed in the 

1950.0 system, x fci , are related to the principal- 

axis (selenodetic) coordinates X ti , by 


x ti ( t) = A T (t)U T (t)X ti 


(4.4.6) 


If p represents the selenodetic radius r, longitude t, 
qj latitude L of an ALSEP , then 


3p 3P 


= a T„T 
~ = 3p 


(4.4.7) 






Since 


h 


*ti 


r cos l cos L 
r sin & cos L 
r sin L 


(4.4.8) 


we have 


ax. . 

— ti 
3r 


iti 


(4.4.9a) 


3 *ti 

30 


-r sin& cos L 
r cos a cos L 
0 


(4.4.9b) 


3X 


•ti 


3L 


-r cosJl sin L 
-r sin H sin L 
r cos L 


(4.4.9c) 


4.4.3 Partials with Respect to Libration 
Parameters 

For physical libration parameters, 

Equation (4.4 i 6) may again be used, since the libra 
tion affects the observable through the lunar rotation 


matrix U(t) 


ax. . ax . 

—X] _ — tx 


a P 


ap 


= A 


T 


au 


T 


= ap -ti 


(4.4.10) 


The elements of are functions of the Euler angles 
ijj, 6, and <j>, which describe the orientation of the 



iiilt-.7;aKij'. -ii. l sssi-nr ti ana. 1 -Ji 
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selenodetic system with respect to the mean-of-date 
ecliptic system (see Section 4.3). Thus, 


9 E T _ *E T a* !E^ ii 

E£> 3p + 36 3p 



3<t> 

3p 


(4.4.11) 


and, from Equation (4.3.9), 


3U 

9 ^ 


-simjjcos4)-cosi|>cosesin(t) cosij>cos4)-sinijjcos0sinij; 0 
sin^sin4>-cosijKX>s0cos<J> -cosi^sincfi— sin4>cos 0cos<() 0 
-cos^sinO -sinijjsine 0 


9 E 

30 


" U 12 U 11 

0 



-u 22 U 21 

0 



~ u 32 U 31 

0 

■ . 



• 

sin^sin0sin<j) 

-cosi^sin0sin(J) 

-oos0sin<J> 

sini|Jsin6cos^ 

-cosij;sin0oos(j) 

-COS0COS(f> 

-sin4cos0 

cosipcos0 

-sin0 


(4.4.12a) 




(4.4.12b) 
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gy -cos^sincf) -s inij>cos 0co s <p -sin^sir^ -hcosipcosdcostp -sin6coscj) 

-~=F' = -c»si|oDS(j)+siniJxx>s0sin({) -sinijxx>s<f)-cosiJraDs9sin<{> sin0sin(|) 


U 21 

U 22 

U 23 


t 

“ U ll 

" U 12 

" U 13 

(4.4.12c) 


0 

0 

0 




lL. i.u 

in which u. . is the element of the l 1 " row and j_ 

* j 

column of U. 

The physical libration parameters affect the Euler 
angles through the libration angles t, a and (I + p) . 
From Equation (4.3.4), 


* = 

+ 

a 



0 = 

I + 

P 


(4.4.13) 

<P = 

d + 

ir 

- + x - a 


so that 





3p 
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9p 



(4.4.14a) 

90 

9p 
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dp 

+ 

9 P 
9p 

(4.4.14b) 

9p 

9t 

3p 
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9a 

3p 

(4.4.14c) 
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The JPL libration tape includes the quantities I 
8 °P 

3 p' and 3^ for the Parameters 3, Y/ C 3Q , C 31 , C 32 , C 33 , 

®31' S 32' S 33' and t ^ le ® initial conditions t q/ p o , Ia Q / 

and Io • We note that the quantity I, which in 


P, 


O' ^O' ~o 

Eckhardt's and other analytic theories is the time-average 
or "mean" value of 0, does not have the same significance 
in the numerically integrated theory. For the numerical 
integration, I was set equal to the mean inclination ob- 
tained by evaluating Eckhardt's theory with the nominal 
parameter values. With £ fixed, the variables p(= 0-1), 

Ict, t, and the partial derivatives of these variables with 
respect to the various parameters were integrated. Thus, 
a partial 3p/3p given on the libration tape is, identically 
80/3p, and 3l/3p = 0. The partials given for (la) and 
(la) are just I(3a/3p) and I(3a/3p), respectively. Note 
also that, by virtue of the definition of I, the mean 
value of p will be non-zero if the parameters are adjusted 
from their nominal values. 

4.4.4 Partials with Respect to Site Coordinates 
and Precession Quantitites 

For parameters which affect the 1950.0 co- 

ax . . 

ordinates of the observing sites, in Equation (4.4.4) 


becomes 


3x . . 

J = 
'3p 


3x . 

- 

3p 


(4.4.15) 


The partials for the geocentric radius, longi- 

tude , and latitude of a site, and for the three angular 
rates which represent corrections to the precession matrix 

are given in detail by Ash [ 6 ] . 

4 .4.5 Partials with Respect to Bias Parameter 

Finally, for the constant "bias " which is added 
to the theoretical value of the observable. 


x 

3P 


(4.4.16) 


for all m > 0. There are no observations processed for 


which m - 0. 
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4.5 Error Analysis 

Any physical process which affects the observations 
and is not properly accounted for in the theoretical 
model of the observable can introduce errors into the 
parameter estimates obtained from analysis of the obser- 
vations. The main sources of these errors may be grouped 
into three general categories: 1) instrumental effects, 

2) propagation medium effects, and 3) deficiencies in the 
dynamical model used to determine the relative motions of 
the ALSEP transmitters and the observing sites. 

4.5.1 Instrumental Effects 

In the class of instrumental effects I in- 
clude those associated with both the transmitters and the 
receiving systems. Of these, the more important effects 
for the observations analyzed in this thesis have been the 
fluctuations of the transmitter frequencies. We saw in 
Section 4.2.3 that under some conditions transmitter fre- 
quency drift could lead to a systematic error in the 
calculation of the theoretical observable exceeding 100° 
of phase. However, in the estimation of parameters I have 
not included any data for which the errors due to frequency 
drift exceeded 5° of phase when averaged over 10 minutes 


or more. 


There is one mode of transmitter frequency fluc- 
tuation which seems to have affected almost all of the 
observations: a quasi-periodic oscillation, with a 

period of 6 to 8 minutes and an amplitude up to 5 Hz*. 

This transmitter-frequency oscillation causes an oscillating 
error in the differenced N-count observable; the ampli- 
tude of the latter oscillation is proportional to the 
propagation-delay difference between the receiving 
stations (see Section 4.2.3). At the beginning or end of 
the tracking period, when the delay difference is largest, 
the amplitude may be as large as 50° of phase, but is 
usually smaller than 25°. Near the mid-time of the ob- 
servation period, the oscillation is imperceptible. Since 
the oscillation is very nearly periodic, it should average 
out in the data even better than errors that are completely 
random. If transmitter-frequency drift is included in the 
calculation of the theoretical observable by means of Equa- 
tion (4.2.22), then the 6- to 8-minute oscillation will be ac 
counted for, provided that the data sampling interval is much 
shorter (~1 minute) than the oscillation period. However, in 
* 

The origin of this oscillation is still being investigated. 
It may be related to cycling of a temperature controller 
for the crystal which determines the transmitter frequency. 
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most of the data processing for this thesis the sampling 
interval was 3 minutes, and Equation (4.2.22) was not used. 

(But see page 83.) 

The dominant sources of error in the receiving systems 
are rate variations of the station clocks*. If the clock 
at one station runs faster than the clock at a second station, 
then the first station will count fewer cycles of the ALSEP 
frequency-difference signal in a given interval of station- 
clock time. The rate of the clock at each station is govern- 
ed by a Cesium-beam primary frequency standard. Tests per- 
formed prior to the Apollo 16 mission [39] indicated that the 
freq uen cies of these standards (or rates of the clocks) may differ 
by at most 1 to 2 parts in 10 12 between different stations. 

If the ALSEP frequency difference is 3 MHz, then a rate off- 
set of 2 parts in 10"^ 2 will cause a variation in the differen 
tial N-count in 6 hours of (3 x 10 Hz)x(2 x 10 ) x 

(2.16 x 10 4 sec) = 0.13 cycles, or 47° of phase. Since the 
effect of a clock rate error is proportional to the trans- 
mitter-frequency difference, it may be possible to identify 
errors stemming from this source by comparing 
the results of simultaneous observations made 


*See also Section 3.2 for a discussion of receiving system errors. 


on the same earth baseline of two ALSEP pairs having 
widely different frequency separations. 

The theoretical observable may also be in 
error as a result of an epoch offset in one 

of the station clocks. The main effect of such an epoch 
error is to cause an error in the calculation of the 
propagation delays. The size of the latter error is 
equal to the epoch error multiplied by the rate of change 
of the difference in delay from two ALSEP transmitters 
to the receiving station. Since the epoch error is un- 
likely ever to exceed 20 usee, and the differenced delay- 

rate is never more than 5 x 10 _9 , the error in the cal- 

-13 

culated delay difference must be less than 10 sec, 
equivalent to less than 0.1° of phase. 

4.5.2 Propagation -Medium Effects 

As discussed in Section 4.2.2, the effect 
of the neutral atmosphere on the propagation delays is 
modeled using monthly average values of the zenith 
delay. The actual atmospheric delay at a station may 
differ from the value calculated from the model by as 
much as 20% as a result of variations in temperature, 
pressure, and water vapor content. We may estimate 
the magnitude of the effect of this error in the ob- 
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servable by using a simple secant law for the dependence 
of the delay on the zenith angle, 0: 


r atm T atm (zenith ) ^ cos ® 


(4.5.1) 


and 


dx 


atm 


de 


— t . , ,sin0/cos 0 

atm (zenith) 


(4.5.2) 


so that the difference in atmospheric delay between 

two ALSEPs separated by A0 in zenith angle is 

2 

Ax . ~ t . , . v (sin0/cos 0)A0. In general 

atm atm (zenith) 

the effect on the differenced N-count observable depends 

on both the orientation of the lunar baseline with respect 

to the horizon at each of the tracking stations, and 

on the elevation of the moon at each station during the 

observing period. However, an estimate of the maximum 

effect can be obtained by calculating the change in 

Ax . at one station as the moon rises (or sets) , as- 
atm 

suming that the other station views the ALSEPs always 
at high elevations, and that the ALSEP baseline appears 
nearly vertical at the former station. The maximum pos- 
sible elevation difference for the most widely separated 
ALSEP pair is 5 x 10 radians. At 10° elevation, the 
value of sin0/cos 0 is~33, so that the amount of dif- 
ferential phase introduced by the atmosphere, if we assume 


-106- 


a zenith delay of 7 ns, is {2.3 x 10 Hz) x (7 x 10 sec) 
x (33) x (5 x 10" 3 ) x (360° /cycle) - 1000°. The effect 
of a 20% error in the model zenith delay is thus about 
200°. However, at 30° elevation, the effect of a 20% 
model error is only 20° of phase. Since the neutral- 
atmospheric delay increases so rapidly near the horizon, 
errors from this source may often be detected from data 
obtained at low elevations. 

Variations of the propagation delay introduced by 
the earth's ionosphere are much more difficult to model 
than those introduced by the neutral atmosphere, because 
the density of the ionosphere fluctuates widely on the 
timescale of hours and days. In order to estimate the 
magnitude of the effect of the ionosphere on the dif- 
ferenced N-count observable, I will use a simple model 
for the elevation-angle dependence of the ionospheric 
delay, together with crude estimates of the time-variabili- 
ty of the ionosphere. I will ignore the detailed de- 
pendences of the ionospheric density on magnetic latitude, 
etc . 

The contribution to the phase delay of a signal 
due to the ionosphere is given by 
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T ion 


^ / N(S)dS 

t path 


(4.5.3) 


where f is the frequency of the signal, N(S) is the 

electron density as a function of position S along 

the propagation path, and k is a constant numerically 

equal to 1.35 x 10” 7 if t ion is in seconds, f in Hz, 

3 

S in meters, and N in electrons/meter . Using a 
modified Chapman model of the ionospheric density 
due to Kazantsev and Nanozov, and assuming a character- 
istic scale height of 350 km, Pfeiffer [26 ] has obtained 
the expression for T ion as a function of the zenith 
angle 8, and of the integrated content at the zenith. 


N„ 


1.35 x 10" 7 N n 


T ion 


H ( 8 ) 


(4.5.4) 


in which the function H(8) takes different forms for high 
and low elevations. With sufficient accuracy for this 
analysis , 

H ( 0 ) =1 + 0.8(sec6 - 1) 0 < 0 * 1 (4.5.5a) 

(35° 5 elevation < 90°) 


H ( 0 ) * 20 - 0.3 


1 < 0 < 1.5 (4.5.5b) 

(5° £ elevation ^ 35°) 


where 8 is measured in radians. 


From Equation (4.5.4) it is evident that there is 
an effect on the observable due both to the frequency 


difference between the ALSEPs and to their geometric 
separation. We will calculate first the geometric 
effect. For nighttime observations, an upper bound on 
N t is 10 17 el/m 2 . Then for two ALSEPs separated by 
A0 in zenith angle 


( Ax ion^ 


1.35 x 10 


10 


dH 


geom 


d9 


■A0 


(4.5.6) 


For elevations between 5° and 35°, dH/d0 ^ 2 


so that an ALSEP pair with an apparent vertical 


-3 

separation of 5 x 10 radians has 


(At. > I ■ I MS X lO™) (2)(5x 10 -3) 
lon (2.3 x 10 y ) 


geom 


= 25 psec 

equivalent to 21 ° of phase. The error introduced into 
the differenced N-count observable at a given time is 
just the difference between the ionospheric effect at 
that time and the effect at the beginning of the track 
(see Section 4.2). However, this difference may be 
expected to be of the order of 100% of the effect 
calculated for low elevation. 
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The frequency effect may be calculated in a similar 
manner. Using Equation (4.5.4), 


dT . 2.7 x 10' 7 N 

= + 3 H<9> 


or 


2.7 x 10" 7 N ( 


(AT ion } 


T 


H (6) Af 


freq 


The maximum frequency separation Af for observations 

with the DDR is 3 MHz. For this frequency separation, 

we find for 9 = 1.2 rad (20° elevation), that 

(Ax. ) « 15 psec, equivalent to 12° of phase, 

ion _ 

freq 

When the moon is more than 60° above the horizon, 

H ( 9 ) ~ 1 and (At. _) * 7 psec. Thus the change 

lon freq 

in delay which affects the differenced N-count, is 
~ 8 ps, equivalent to 7° of phase. 

We may expect that ionospheric effects will be 

more serious for observations of the ALSEP 12-ALSEP 17 

-3 

pair, which has a large geometric separation (5 x 10 rad) 
as well as the maximum possible frequency separation 
(3 MHz) , than for other pairs. For observations of the 
12-17 pair, the total error introduced by the nighttime 
ionosphere may be as great as 30° of phase. However, 


LI 


in the daytime, the total ionospheric electron content 
may. rise to 10 times the nighttime value, with most of 
the increase taking place over a period of 4 hours fol- 
lowing sunrise. For daytime observations especially, 
and for the ALSEP pairs which have the greatest geometri 
cal and frequency separations, the ionospheric effect on 
the observable may produce errors of the order of 
hundreds of degrees of phase. 

Any model for the ionosphere which is not based on 
direct measurements of the electron content on an ap- 
proximately hourly basis seems likely to be in error 
by 50% or more on most days. Consequently, I have 
chosen not to use the relatively crude model available 
in PEP in the processing of the data for this thesis. 
However, in Chapters V and VI I discuss the probable 
effect on the parameter estimates of errors introduced 
by the ionosphere and how eventual reprocessing of the 
data with an improved model might reduce the effect of 
these errors. 

4.5.3 Effect of Errors in the Dynamical 
Model of the Earth-Moon System 

There are two obvious requirements 

for the proper use of a dynamical model in fitting 

a set of observations: first, the model itself must 
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accurately represent the dynamics, and second^ all the 
parameters in the model which affect the observations must 
be adjusted simultaneously in the least-^squares estimation 
procedure. For ALSEP differential VLB! observations, these 
two requirements apply primarily to the model of the moon's 
physical libration, but also to the models for the moon's 
orbital motion and the earth's rotation. 

In considering possible errors in PEP's model of the 
dynamics of the earth-moon system, I have relied heavily on 
the results of fits I have made to the lunar laser ranging 
data. In these fits I adjusted the moon's 6 orbital elements 
the mass of the earth-moon system, up to 15 parameters affecting 
the physical libration, and 3 corrections to the adopted 
precession matrix. These fits are useful because the laser 
ranging data are more sensitive than the differential VLBI 
data to most of •'.he dynamics, especially the orbital motion 
of the moon and the earth's rotation. An indication of the 
size of the errors in the latter models is afforded by the 
fact thcit I have been able to fit the laser data with rms 
residuals of 6Q cm in (one-way) range. It is possible, of 
course, that errors much greater than 60 cm exist in the 
model but are being masked in the fitting process. Since 
the laser observations measure- range , and are relatively 
insensitive to rotational motions of the moon about the 
line of sight, there may be particularly large - orientation 
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errors in libration models which are derived from fits to 
laser data alone, 

• The ALSEP differential VLBI data are more sensitive to 
the moon's libration than to the orbital motion. The model 
which I used for the libration was the JPL LLB— 5 numerical 
integration. Williams et al . [41] have tested for errors 
in the integration by comparing it with a semi-analytic 
libration model developed independently by Eckhardt [15,16]. 
When a 7-year integration of the libration angles p , Ict , 
and x is fit by least-squares to the Eckhardt series "300" 
libration model, the rms difference between the two models 
for these angles is less than 0','2, equivalent to approximately 
1.6 m at the lunar surface. 
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CHAPTER V 
RESULTS 


5 . 1 Introduction 

In this chapter I present the results of solutions 
for the libration parameters and the ALSEP coordinates. 

AH of these solutions used combinations of data from 
differential VLBI observations of the five ALSEP trans- 
mitters, together with data from laser ranging observa- 
tions of the Apollo 11, 14, and 15 Lunar Ranging Retro- 
reflectors (LRRRs) . The laser ranging observations were 
performed at the Me Donald Observatory under the auspices 
of the NASA Lunar Laser Ranging Experiment (LURE) [ 7 ] . 

The data obtained from the laser observations are in the 
form of 1194 normal points constructed at the University 
of Texas from ranges obtained between January 1970 and 
June 1974 [ 1 , 33 , 34] . 

The VLBI data were selected from observations made be 
tween March 1973 and June 1974 at the six STDN stations 
which are equipped with DDRs (see Figure 1.5)*. These 
observations involved all five ALSEP transmitters, in 
8 pair-wise combinations.** 


*One of these 6 stations, the Network Test and Training Faci- 
lity at Greenbelt, Md., was actually used for ALSEP DDR observations 
on only 2 days: 7 and 8 March 1973. 

** Of the 10 pairs which can be formed, from the 5 ALSEPs, 2 P^irs 
(AL14-AL17; AL14-AL16) cannot be observed byfheDDR because e 
transmitter-frequency separations are more than 3.3 MH. (see sec 
tion 3.2) . 
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Since the VLBI observations are sensitive primarily 
to the relative positions of |the different ALSEPs, I have 
used an independent photogrammetric determination of the 
coordinates of the Apollo 15 ALSEP with respect to the 
nearby LRRR, in order to tie the system Jf ALSEP co- 
ordinates to that of the laser retroref lectors. The 
photogrammetric determination was done by Schimerman et al. 
of the Defense Mapping Agency Aerospace Center, using 
the photographs of the landing site taken by the astro- 
nauts. The azimuth and the distance of the Apollo 15 
ALSEP from the LRRR were obtained with uncertainties of 
about 1 meter [30 ] . 

In Appendix A. 2, I have listed all of the ALSEP ob- 
servations made through June 1974 and have designated 
those which were used in the solutions. The observations 
are grouped into series , each consisting of a sequence 
of observations of a pair of ALSEPs by a pair of track- 
ing stations during which the cycle-count at both sta- 
tions was uninterrupted. Since each STDN station is 
equipped with two DDRs, and since more than two stations 
often observe at the same time, several observation 
series were often obtained on a single day. As discussed 
in Sections 2.2, 4.1, and 4.2, each observation series 
has associated with it an adjustable bias parameter 
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which accounts for the error in the initial observation 
in the series. At present, PEP can accommodate only 
130 biases, and hence only 130 observation series, in 
a single solution. For this reason, as well as to pro- 
vide a check of my results, I divided the ALSEP data into 

two groups. 

The first group consists of 10,602 observations 
(100 observation series) of ALSEP 15 with respect to 
other ALSEPs. These observations have been combined with 
the 1194 laser normal points to form the data set which 
I will designate A (DS A) . The second group consists of 
the remaining observations (10,123) and includes obser- 
vations of all pairs which can be formed from the 
other four ALSEPs. This group, combined with the laser 
normal points, comprise Data Set B (DS B) . The particular 
division of the data which I have chosen has several 
advantages. First, the inclusion of all observations 
of ALSEP 15 in DS A allows a determination of the co- 
ordinates of ALSEPs 12, 14, 16, and 17 with respect to 
the Apollo 15 retroref lector . Second, because ALSEP 15 
is absent from DS B the geometry of the VLBI observations 
in this data set is significantly different from that of 
DS A (see Figure 5.1) , so that a meaningful consistency 






of the ALSEP co- 


check is possible. The origin 
ordinate system for the DS B solutions was defined ar- 
bitrarily by fixing the coordinates of ALSEP 16 at values 
which had been determined by preliminary solutions and 
which are probably accurate to within 30 m. The dis- 
advantage of this grouping of the data is that the para- 
meter estimates obtained using DS B may have errors in- 
troduced by the errors in the assumed ALSEP 16 coordinates 
For this reason, a third data set was formed, DS C, 
consisting of DS B plus 1569 observations of ALSEP 16 with 
respect to ALSEP 15 taken from DS A. A summary of the 
observations contained in each of the three data sets 
is given in Table 5.1. 

With each of these three data sets I performed a 
series of solutions. In every solution I estimated the 
6 elements of the lunar orbit, the mass of the earth- 
moon system; three angular rates representing corrections 
to the earth's precession matrix; 3 coordinates each for 
the McDonald Observatory, the STDN stations*, the Apollo 


A priori information was included for each of the STDN 
station coordinates as shown in Table 5.2. 
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Table 5 . 1 

Differential VLB! Observations Contained in the 

Data Sets* 


ALSEP Pair 

No. Obs. 

No . Obs . 

No. Diff 

. Period 




Series 

Tracking 

Days 

Covered 


Data Set A 







AL 12-AL 15 

3418 

34 * i 

24 

Mar 

73-Jun 

74 

AL 14-AL 15 

2153 

20 

14 

Mar 

73-Jun 

74 

Al 16-AL 15 

2719 

29 

18 

Jul 

7 3 -May 

74 

AL 17-A1 15 

2312 

17 

11 

Apr 

73-Jun 

74 

Totals : 

10,602 

100 

49 

Mar 

73-Jun 

74 

Data Set B 







AL 12-AL IS. 

3509 

36 

27 

Jul 

73-Jun 

74 

Al 12-AL 16 

2281 

23 

18 

Aug 

73-Jun 

74 

AL 12-AL 17 

2781! 

18 

14 

Jul 

7 3 -Feb 

74 

AL 16-AL 17 

1552 

8 

7 

Nov 

7 3 -May 

74 

Totals : 

10,123 

85 

48 

Jul 

73-Jun 

74 

Data Set C 







From Data 
Set B 

10,123 

85 

48 

Jul 

73-Jun 

74 

AL 16-A1 15 
from Data 
Set A 

1569 

17 

12 

Sep 

7 3 -May 

74 

Totals : 

11,692 

102 

53 

Jul 

73-Jun 

74 


* 

Each Data Set also contains 1194 laser ranging normal 
points covering the period Jan 1970 - June 1974. 




Table 5.2 


Nominal STDN Station Coordinates* 


Station 


Ascension Is. (AON) 

Madrid, Spain (MAD) 

Merritt Is., Fla. 

(MIL) 

Merritt Is. , Fla. - 
System 2 (MILW) 

Greenbelt, Md. (ETC) 

Corpus Christi, Tex. 
(TEX) 

Goldstone, Calif. 
(GDS) 


Geodetic Coordinates ** 

Height Long itude Geodetic Lat. 
(m) (deg) (deg.) 


527.0 

14.327058 

-7.954794 

763.0 

4.168453 

40.454892 

-54.0 

80.693469 

28.508272 

-54.0 

8Q. 693469 

28.507753 

o 

• 

oo 

CN 

1 

76.842914 

38.998481 

-39.0 

97.378622 

27,653750 

921.0 

116.873597 

35.341594 


♦The a priori standard deviations for Ascensioij and Madrid 
were 15 m in geocentric radius, and 1.5 x 10*"' deg 
longitude and latitude; for all other stations, 10 m xn 
radius, and 1,0 x 10 -4 deg in longitude and latitude were 
used. 

♦♦Referred to an ellipsoid with equatorial radius 

6378.166 km and degree of flattening 1/298,3. See [6], 




retroref lectors , and the ALSEP transmitters; a bias 
parameter for each VLBI observation series; the lunar 
moment-of-inertia parameters $ and y ; certain third- 
degree-harmonic coefficients in the lunar gravitational 
potential; and the 6 initial conditions of the lunar 
libration. Different combinations of the third-degree- 
harmonic coefficients were solved for in the different 
solutions. 
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Figure 5.2 shows examples of the postfit residuals 
obtained from the solutions. The best fits of the VLBI 
data were obtained for observations of ALSEPs 12 and 14 
[graph (a)]. These ALSEPs are separated by only 180 km 
on the moon, so that the observations are relatively 
insensitive to the physical libration and also to the 
effects of the earth's atmosphere and ionosphere. Graph 
(b) shows one of the worst^f itting series (186D-) and tvro 
others (186F and 186E) whose r,m,s. residuals were also 
larger than the average for all of the VLBI data (16°). 

Since a bias parameter was adjusted for each series, errors 
in the theoretical model tend to show up primarily as a 
slope in the postfit residuals. The slope in the residuals 
for series 186D (~40°/hour) is approximately equivalent to 
a 10 meter error in ALSEP relative position. 

The laser residuals shown in the bottom graph of 
Figure 5.2 are typical for these solutions. The r.m.s, 
of the residuals for this 8-day span is 2,7 ns, compared 
with 3.5 ns for all of the laser data. The random noise 
level for the normal points in this span, based on the 
scatter in individual photon returns, is ~0.8 ns. The 
residuals shown in the graph illustrate the results of three 
different types of model error. The 1-2 ns trend which is 
common to all three reflectors results from lunar orbit 
errors which affect the earth-^moon center-^to^-center 


.(a) VLB I POSTFIT RESIDUALS 28 SEPT 1973 

© Series 271C ACN-MAD ALSEPs 12 & 14 13 rms) 
m a Seri es 271D ACN-MIL ALSEPs 12 & 14 ( 7 rms) 
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\ . 


VLB I POSTFIT RESIDUALS 5-6 JULY 1973 
© Series 186n MAD -AC N ALSEPS 15 & 16 (27° rms) 

* Series 186F MAD-MIL ALSEPs 15 & 6 <iq° rms) 

+ Series 186E MIL-GDS ALSEPs 15 & 16 (19 Arms) 

▲ 
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(c) LASER RANGING POSTFIT RESIDUALS 4-11 NOV 1973 


A Apollo 11 
+ Apollo 14 
O Apollo 15 



^ ST 1794 1795 1796 1797 1799 1799 

Julian Day Relative to 2440000 

5 2 Typical postf it residuals from solutions combining both 
VLB I and laser ranging data. The upper . wo grap 
examples of the best (a) and worst (b) fits of VLBI 
series. The overall rms of the VLBI residuals was 16 . 
The bottom graph (o) is typical of the laser ranging 
residuals throughout. (See text.) 
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distance. Differences in the residuals between different 
reflectors at the same time, also about 2 ns, can be 
attributed mainly to errors in the physical libration 
model. Finally, a spread in the residuals obtained from 
a single reflector over one day may result from an error 
in the model for earth rotation. On two of the days 
shown, the spread exceeds 5 ns. 

5.2 Libration Parameters 

The values of the libration parameters used by 
Williams et al , at JPL to generate the numerically- 
integrated libration model, LLB-5, were derived par- 
tially from fits to laser ranging data and partially 
from the results of solutions performed by other groups 
using Doppler tracking data from the Lunar Orbiter space- 
craft [43]. In particular, the values used for J 2f C 31 , C 33 
and s 3 ^ were calculated from an average of the Lunar 
Orbiter results of Michael and Blackshear [22] , Liu and 
Laing [21], and Sjogren [35]; 3, Y, and the remaining 
third degree harmonics were determined from solutions 
performed at JPL using essentially the same set of laser 
normal points as I have used in my solutions*. Prior to 

*The time period covered by the normal points used at JPL 
and at MIT is identical. Minor differences exist between 
the two data sets as a result of different editing criteria 
and the fact that the JPL data set includes some ranges to 
Lunakhod II [43] . 
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performing solutions with the ALSEP differential VLBI data, 

I used the laser ranging data alone in an effort to duplicate 
the solutions performed at JPL, The estimates obtained for 
libration parameters are given in columns 3 and 4 of Table 
5.3. With a nondegenerate set of harmonic coefficients 
ydLimdtsd (column 3 of the table , PS— 1) , most q f m^ results 
are consistent with the nominal values of LLB-5. 

The first solutions I performed with Data Set A indicated 
that the addition of the VLBI data to the laser data improves 
substantially the formal errors obtained for C 31 and C 33 , 
which were poorly determined with laser data alone, (Results 
with laser data alone are shown for two parameter sets, one 
[PS-1] in Which C 31 , C 33 , and S 31 were fixed and one [PS-2] 
in which all harmonics were estimated.). The improvement in 
the formal error for C 31 is by a factor of 100 , and for C 33 , a 
factor of 25. These improvements are also reflected in the 
parameter correlations. With the laser data alone, the 
correlation coefficient between C 33 and y is r-, 99.7. The 
addition of the VLBI data reduces the correlation to -0.92. 
With laser data alone, C 31 has a correlation of 0.97 with 3; 
with VLBI data added, the correlation coefficient is 0,003, 

The formal errors shown in Table 5,3 do not necessarily 
reflect the true uncertainties in the parameter estimates. 

The postfit residuals for both laser and ALSEP data are not 
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Table 5.3 

Comparison of Libration Parameter Values and Formal Standard 
Errors* Obtained with Laser and VLBI Data Sets 

(Units of 1(T 6 ) 


Parameter 

Nominal 

Value 

Laser Data 

Alone 

Laser + 

VLBI Data 


(LLB-5) 

PS-lt 

PS-2tt 

(DS A) 

PS-lf PS-2 ft 

e 

631.26 

631. 28±0 . 01 

629±1 

631. 266±0 , 003 

631 . 261±0 ,003 

Y 

227.37 

227 . 39*0.03 

224±1 

227 , 42±0 . 01 

227 , 34±0 . 04 

J 3 

10.44 

5±2 

7±2 

5. 6±0. 8 

1±1 

C 31 

28. 6t+f 

28.6 

21±20 

28.6 

28 , 3±0 . 2 

C 32 

4.82 

4 , 75±0 . 04 

4 .70±0 , 05 

4 . 83±0 . 02 

4 . 70±0 . 03 

C 33 

2.7 

2.7 

9±2 

2.7 

2 . 82±0 .08 

S 31 

8.8 

8.8 

12±7 

8.8 

-8±4 

S 32 

1.71 

1. 66±0, 06 

1. 45±0 . 08 

1 . 7 2 ± 0. , 0 3 

1. 79±0 . 04 

S 33 

-1.14 

0 . 0±0 . 2 

0 . 0±0 . 2 

-0 , 2±0, 1 

-0 . 4±0 . 1 


♦Formal standard errors (based on the r.m.s. values of the postfit 
residuals) are given. here for the purpose of comparing 
the sensitivities of the two data sets to the parameters. The actual 
uncertainties of the parameter values are typically much larger 
(see Table 5.5), 

tParameter Set 1 included 6 initial conditions of the lunar orbit, 
mass of earth + moon, 3 adjustments to the earth's precession matrix, 
3 coordinates each of McDonald Observatory and the Apollo 11, 14, and 
15 LRRRs , 6, Y, J 3 , C 32 , S 32 , S 33 and 6 initial conditions of the 
libration. When VLBI data were used the parameter set also included 
3 coordinates each of the 6 tracking stations (with 10—15 m a priori 
a's) and 4 ALSEP transmitters (AL 15 fixed to AP 15 LRRR) , and a 
bias parameter for each of 100 observation series. 

f tParameter Set 2 is the same as PS-1 with the addition of S 3 i, C 31 , 
and C 33 . 

tttUnderlined values were not adjusted to fit either laser or VLBI 
data but were fixed from an average of Lunar Orbiter results (see 
Table 5.5 ) . 
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j • random. In order to obtain a fairer estimate of the uncer- 

| tainties in the parameter values , I performed a series of 

| 

\ ' solutions in which different combinations of parameters 

| were estimated with different data sets. For each parameter 

| set a solution was first made with laser data alone, then 

| 

j w ith Data Sets A, b, and C, each of which includes combma- 

i tions of laser and VLBI data. Table 5.4 s hows a comparison 

1 

jj of the results obtained with Data Sets A, B, and C when all 

f: 

| the libration parameters were adjusted. The agreement among 

| ' the three solutions is poorest (inconsistency on the order of 

| 20 times the formal standard error) for the three parameters 

I ( Y, C 31 , and C 33 ) which are most highly correlated, possibly 

reflecting a larger influence on these parameters of systematic 

. . ' : 

i errors . 

;i 

I My best estimates of the values of the libratio.n param - 

eters and their uncertainties are summarized in Table 5.5. 

For comparison, I have also included in the table three deter- 
minations of the third-degree harmonic coefficients from Lunar 
Orbiter Doppler tracking data, and the determinations from 
lunar laser ranging data reported recently by J.G. Williams 
of JPL. The VLBI results listed are averages of the values 
obtained from the DS A and DS C solutions.* The uncertainties 


*The DS C solution is believed to be more reliable than the 

DS B solution because, in the former, VLBI observations of 
the Apollo 15 ALSEP were included in order to relate the 
ALSEP positions to those of the laser retroref lectors . 
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Table 5.4 

Comparison of Solutions for Libration Parameters 
with Data Sets A, B, and C 
(units of 10' 6 ) 


Parameter 

DS A 

++ 
DS B 

DSC 

3 

631 . 261±0 . 003 + 

631. 27310. 003 

631,277+0.003 

Y 

227. 34 ±0,04 

223.09 +0.04 

228. 0510. 04 

J 3 

1±1 

6 ±1 

511 

C 31 

28.3±0.2 

24.4+0,2 

24.510.2 

C 32 

4.70±0.02 

4.8410.02 

4.7710,02 

C 33 

2 . 810.1 

1,310.1 

1, 3+0.1 

S 31 

-8 ±4 

26+3 

5 14 

s 32 

1.7910.04 

1.7910.03 

1.7610.03 

S 33 

-0.410.1 

-1. 010. 1 

-0,210.1 

+ Formal standard error 
++ ALSEP 16 coordinates fixed: 

Radius = 1737. 

447 km, 


Long. = 15.43625°, Lat. = -8.95541° 
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Table 5.5 

Estimates of Lunar Gravity and Libration Parameters 

(units of 10 


Lunar Orbiter Doppler Laser Ranging 


Parameter 

Michael 
& Black 
Shear 
[22] 

Liu & 
- Laing 
[21] 

Sjogren 

[35] 

Wi 1 ] i ams 
[42]* 

+VLBI* 

3 

— 

■ 

— 

631.2610.3 

631.27+0.03 

Y 

- — ■ 

r 


227.3710.6 

227.7+0.7 

J 2 

2071, 

1996±2 

204813 

2038. 22 + 

2038. 22 + 

C 22 

22.4 

23. 6+5.3 

22.110.5 

22 . 4 ++ 

22 , 4 ++ 

J 3 

-28.4 

5 . 9±2 . 9 

10.719 

10.44+4 

3i20 

C 31 

24.1 

30 . 0±3 . 7 

31.615 

28. 6 + 

2614 

C 32 

7.6 

4.712.8 

5.511 

4.8210.15 

4.7+0.2 

C 33 

1.4 

4.812.2 

1.810.3 

2.7 

2 +2 

S 31 

20.8 

1.413,2 

4.3112 

8.8 

-1 + 30 

S 32 

2.27 

0.611.7 

2.7 + 1 

1.7110.15 

1.8 + 0.3 

S 33 

-0.3 - 

■0.9 + 1.3 

-1.010.2 

-1.14+0.7 

-0.311.0 


♦Estimated uncertainties are provisional 

+ Underlined values for laser and VLBI solutions were fixed 
from an average of the Lunar Orbiter determinations 

++ Derived from 3/ Y* and J 
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assigned to most of my results are a factor of 20 larger than 
the corresponding formal standard errors. This factor is 
based on the observed disagreements among the DS A, B, and C 
solutions in Table 5.4. For C-^r ^31' ^32' an< ^ ^33 * have 
used 10a rather than 20a uncertainties since these coefficients 
are determined almost entirely by the laser data (see Table 
5.3) , which I believe to be less corrupted than the VLBI data 
by systematic errors. 
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5 . 3 ALSEP Coordinates 

The coordinates of the ALSEP transmitters were adjusted 
in all solutions which included VLBI data. The results to be 
presented in this section were obtained from the solutions with 
PS-2, for which all of the libration parameters were estimated. 
Since the third-degree harmonics produce offsets of the mean 
directions of the moon's principal axes from the mean direction 
of the earth, the selenodetic (principal-axis) coordinate systems 
defined by solutions obtained using different data sets tend to 
differ. In an effort to define a consistent selenodetic co- 
ordinate system for the comparison of results from solutions with 
Data Sets A, B, and C, in all of these solutions I fixed the co- 
ordinates of the Apollo 15 LRRR at the values determined from the 
laser-only solution. Fixing the Apollo 15 LRRR coordinates in 
this way reduces substantially the correlations among the retro- 
reflector and the ALSEP coordinates, but has only a negligible 
effect on the r.m.s. of the post-fit residuals. 

Since the primary tie between the ALSEPs and the laser 
retroref lectors is at the Apollo 15 site (see Section 5.1), 

Data Set A (which includes all of the data from observations 
of ALSEP 15) was used in a solution to determine the 
coordinates of the other ALSEPs with respect to the Apollo 15 
LRRR. The results of this solution are given in Table 5.6. 

These results were then checked by performing additional 
solutions for the relative coordinates of the Apollo 12, 14, 

16, and 17 ALSEPs using Data Sets B and C. For solutions with 
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Table 5.6 

3 

Relative Coordinates (ALSEP-Apollo 15 LR ) Using 
Adjusted Libration Model* 


ALSEP 

Mi 

(km) 

ax 2 

(km) 

AX 3 

(km) 

ARadius 

(km) 

ALong. 

(deg) 

ALat . 
(deg. ) 

12 

35.287 

-787.824 

-855.543 

0.503 - 

■27.052 9 6 

-29.14417 

14 

97.589 

-619.135 

-874.730 

0.860 - 

■21.10 591 

-29.77 784 

15** 

- 0.012 

0.038 

0.020 

0 

0.00143 

0.00074 

16 

99.611 

359.859 

-494.678 

1.9 44 

11.86812 

-35.10912 

17 

-155.057 

734.404 

-165.681 

-0.695 

27.13959 

- 5.94397 


'From solution with Data Set A. The Apollo 15 LR 3 selenodetic 
coordinates were fixed at the following values^ 40 

Radius = 1735.5087 km, Long. = 3.56822 , -at. 26.15 

t- P 2 - 7 - 

® 31 = -2oi'8i, S ia o = 1 26?’i2, 33 T 0 = -5V492/day, 0 O = -l'.'938/day, 


r o 

la. 


= -5'.'04 4/day. 


**Fixed in the solution. These values are from the photo- 
grammetric determination by Schimerman et a l. [30]. 
uncertainty in each coordinate is 1 m. 
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DS B, which contains no data from observations of ALSEP 15, 
the ALSEP coordinate system was fixed relative to the retro- 

reflectors by fixing the 3 coordinates of ALSEP 16, as described 

$ 

in Section 5.1. The ALSEP relative coordinates determined 
from the three solutions, with Data Sets A, B, and C, 
are in Table 5.7. For the short (180 km) ALSEP 12 - 
ALSEP 14 lunar baseline, the agreement among the 3 deter- 
minations is within 3 meters in all 3 coordinates. 

For the other baselines, the results of the 3 solutions 
disagree by up to 20 m in the radial (X-^) coordinate and by up 
to 8 m in the 2 transverse coordinates (X 2 and X 3 ) . 

An important external check of the VLBI results is 
also available. The relative coordinates of the Apollo 14^ 

ALSEP and LRRR determined from the landing-site photo- 
graphs may be compared with the coordinate differences 
determined from the earth-based observations. The co- 
ordinates of the Apollo 14 LRRR with respect to the 
Apollo 15 LRRR have been determined by laser ranging, 
and the coordinates of ALSEP 14 with respect to ALSEP 15 
have been determined via differential VLBI*. If we 


These two determinations are not entirely independent 
inasmuch as both were based on the use of the same 
libration model. 
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Table 5 . 7 

Comparison of Solutions for ALSEP Relative Coordinates 
with Adjusted Libration Model 


ALSEPS 

Data 

Set 

AX, 

(km) 

Dif f . 

from 

Soln. 

..A 

ax 9 

(km) 

Dif f . 
from 
Soln. 
A 

AX, 

(km) 

Diff . 
from 
Soln. 
A 




(m) 


(m) 


(m) 

12-14 

A 

- 62.304 

— - 

-168.688 

— 

19.190 

— 


B 

- 62.306 

- 2 

-168.688 

0 

19.192 

2 


C 

- 62.304 

0 

-168.689 

-1 

19.192 

2 

16-14 

A 

2.013 

— 

978.995 

— 

-160.704 

— 


B 

1.994 

-19 

979.000 

5 

-160.708 

~4 


C 

1.994 

-19 

979.002 

7 

-160.708 

-4 

17-14 

A 

-252.619 

- — 

1353.535 

— 

709.061 

— 


B 

-252.618 

1 

1353.542 

7 

709.054 

-7 


C 

-252.611 

8 

1353.542 

7 

709.053 

-8 

17-16 

A 

-254.632 

— 

374.541 

* 

869.765 

— 


B 

-254.612 

20 

374.542 

1 

869.762 

—3 


C 

-254.614 

18 

374.542 

1 

869.761 

-4 
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assume that the relative coordinates of the ALSEP trans- 
mitter and LRRR at the Apollo 15 site are well determined 
(uncertainties less than 1 m) from photogrammetric 
analysis, then a comparison of the relative coordinates 
of the LRRR and ALSEP at the Apollo 14 site determined 
by laser ranging and VLBI with the relative coordinates 
determined from the landing site photographs constitutes 
an excellent check on the accuracy of both earth-based 
observation types. The results of this comparison are 
given in Table 5.8. The disparity between the coordinates 
determined by earth-based and moon-based methods is only 
8 m and may in fact include errors at the level of 5 m 
in the as yet preliminary photogrammetric determination 
at the Apollo 14 site. 
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Table 5 . 8 


Comparison of Earth-based and Moon-based Determinations 
of the Relative Coordinates of the ALSEP Transmitter with 
Respect to the LRRR at the Apollo 14 Landing Site 


Method 

Difference (ALSEP-LRRR) 

Distance 

Longitude 

Latitude Radius 



Laser Ranging + 
VLB I* 

0°. 00101 

0°. 00047 7 m 

34 

m 

Apollo 14 
Landing Site 
Photographs* * 

0°. 00132 

0°. 00023 

41 

m 

Difference 
(Laser/VLBI 
- Photo. ) 

-0°. 00031 
(-8 m) 

0°. 00024 

(7 m) 

-7 

m 


*from solution with Data Set A. In this solution the 
relative coordinates of ALSEP 15 with respect to the 
Aoollo 15 LRRR were taken from the photogramme trie 
determination by Schimerman et al . t 3 0 1 • See 
Table 5.6. 

**from preliminary results with uncertainties of the 
order of 5 m [ 8] . 
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5.4 Discussion 

The estimates of uncertainty which have been derived in 
the preceding sections from comparisons of different solutions 
with different data sets and parameter sets may not be 
completely reliable. Such comparisons do not necessarily 
uncover every kind of systematic error that might be present. 
One reason for caution is the presence of obvious systematic 
trends in both the VLBI and laser residuals (see Figure 5.2). 
My results may also have been affected significantly by 
errors in the partial derivatives on the libration tape. 

These partials were generated from earlier numerical integra- 
tions, which were based on different values of the parameters. 
Consequently, the errors in some of the partials may be of the 
order of 1% [43], Errors in the partial derivatives tend to 
affect most seriously the estimates of parameters which are 
highly correlated with other parameters. However, because I 
was unable to re-integrate the libration using the adjusted 
parameter values, I cannot say just how serious the effects 
of partial-derivative errors were. 

Additional questions have arisen as a result of comparing 
my laser-only solutions with similar solutions performed at 
JPL by J.G. Williams and W.S. Sinclair [43]. With the 
libration model fixed (i.e., LLB— 5, unadjusted) , I obtain 
retroref lector coordinates which differ from those obtained 
at JPL by 10-15 m. The coordinates estimated for the 
McDonald Observatory also differ, by 0.°0013 (12 m) in 
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longitude, 0.1 m in cylindrical radius, and 2 m in the z 
coordinate. The discrepancy for McDonald may be due to 
differences between the JPL and MXT planetary ephemerides , 
which effectively determine the orientations of our respective 
coordinate systems, I attempted to absorb possible orientation 
errors by adjusting three angular-rate corrections to the 
precession matrix (see Section 4.2.2 and also Ash [6] , p. 60) . 

The estimated rates were of the order of l K /century about all 
three axes. The STDN station coordinates were also adjusted 
(with 10-15 m a priori a's), partly to account for any errors 
in orientation of the network with respect to McDonald, 

However, the adjustments to these coordinates were sometimes excessive 
(up to 200 m in the radii and latitudes of Ascension and Madrid) 
and showed little consistency between the DSA and DSB solutions. 
This inconsistency may indicate that the STDN station coordinates 
were more sensitive than other parameters to unmodeled effects 
of the ionosphere and atmosphere. 

I expect that many of the problems considered in this 
section will be better understood, if not eliminated, as the 
relevant models are improved and additional data become available. 
Work on the models is continuing at MXT and elsewhere. The 
observations also continue, with certain improvements as 
discussed in Chapter 6. 
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CHAPTER VI 
CONCLUSION 

6 . 1 Summary 

Earth-based radio observations have been used for the first 
time to determine the positions of transmitters on the surface 
of the moon with uncertainties less than 30 m, and to determine 
the moon's. physical libration. In these observations a new 
phase-locked differential tracking technique was employed, and 
an interferometric observable was formed by subtracting 
simultaneous measurements at widely separated ground stations. 

The differential interferometric data were combined with 
lunar laser ranging data to estimate parameters of the moon's 
physical libration and to determine the relative seleno centric 
coordinates of the 5 Apollo Lunar Surface Experiments Package 
(ALSEP) transmitters. The parameter estimation was carried out 
within the framework of the MIT Planetary Ephemeris Program, 
using a numerically-integrated physical libration model developed 
at the Jet Propulsion Laboratory. The estimate obtained for 
the lunar moment-of-inertia ratio 0 is believed to be more 
accurate than any previous determination. The uncertainties 
of the estimates obtained for 2 of the third-degree harmonic 
coefficients of the moon's gravitational potential are smaller 
than those obtained using laser ranging data alone and are 
comparable to the uncertainties obtained by other investigators 
using Lunar Orbiter Doppler data. The values obtained for all 
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of the harmonics are consistent with the Lunar Orbiter results. 
The relative coordinates of the ALSEP transmitters were 
determined with uncertainties of about 30 m in the radial (X^) 
coordinates and about 10 m in the two orthogonal coordinates 
(X 2 and X^) . The dominant sources of error in all of the 
parameter estimates may have been unmodeled effects of the 
propagation medium. 

6.2 Suggestions for Further Study 

Future work should be directed toward both the improvement 
of dynamical models and the reduction of observational errors. 
The modifications necessary to enable PEP to integrate directly 
the equations of motion for the rotation of the moon are 
already in progress. The completion of these modifications 
will provide an independent check of the JPL integration. 
Further, it will be possible to re-integrate the libration 
using adjusted parameter values, in order to minimize the 
effects of errors in the partial derivatives (see Section 5.4), 
The moon's elasticity (and inelasticity) should be accounted 
for in the equations of motion for the rotation. 

The three sources of error believed to be most serious 
for the VLBI observations analyzed in this thesis - the 
ionosphere, the neutral atmosphere, and variations in the rates 
of the station clocks - may be significantly reduced for the 
observations now being made with improved equipment, A new 
device, called the Differential Doppler Processor (DDPV allows 


the use of multiple DDRs at a station, so that all 5 ALSEPs 
may be observed simultaneously. The additional information 
available at each sample time may be used to construct a 
new observable, a linear combination of the 4 differenced 
N-counts, which has minimum sensitivity to certain carefully 
defined effects; for example, effects which are dependent on 
elevation above the horizon at a station (atmosphere and 
ionosphere) and effects which are dependent on frequency 
(ionosphere and station clock rate variations) , 

For the existing VLBI data, additional study is needed 
to determine whether reprocessing with improved modeling of 
the ionosphere would yield significantly improved estimates 
of the libration parameters and the ALSEP coordinates, 

6.3 Applications 

The results of this work have immediate or near-term 
application in at least four areas of astronomical research; 

1) lunar mapping, 2) studies of the moon's gravitational field, 

3) astrometry, and 4) radio tracking of future lunar and 
planetary spacecraft. 

The ALSEP relative positions presented in Table 5,6 are of 
sufficient accuracy (10-30 m) to be of immediate use in providing 
calibration points for combining orbital metric photography from 
three different Apollo missions (15, 16, and 17) [31], Errors 

in the models for the Apollo spacecraft orbits have limited 
the accuracy with which photographic features could be located 
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in a selenocentric coordinate system to the level of several 
hundred meters, an order of magnitude larger than the resolution 
limit of the photographs. The inclusion of the accurate 
coordinates for the ALSEPs in the photogrammetrrc solutions 
should help to bring the uncertainties closer to the resolution 

limit. 

The values estimated in this thesis for 6, Y, and the 
third-degree gravitational harmonic coefficients (Table 5.5) 
represent only a slight improvement over the results from lunar- 
orbiter tracking and laser ranging. However, the sensitivity 
of the VLBI observations to 6,C 31 , and C^, indicated by 
formal errors, suggests that the values of these parameters 
will be improved further as more data become available. 

The results of this experiment may be important in 
future uses of the moon for astrometric and geodetic purposes. 
With lunar laser ranging and with the photogrammetric tie 
between the Apollo 15 retroref lector and ALSEP transmitter, 
the VLBI results determine the positions of the ALSEPs 
with respect to the moon’s center of mass. The uncertainty 
in the selenodetic positions, set mainly by the uncertainty 
in parameters in the libration model, is approximately JO m, 
corresponding to an angular uncertainty of 0V015 of geocentric 
arc. With this tie to the lunar center of mass, the ALSEPs 
may be used in differential VLBI observations between the moon 
and extragalactic radio sources in order to determine the 
orientation of the lunar orbit with respect to an inertial 
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reference frame [11] . A program of such observations is presently 
being carried out by M.A. Slade and R.A. Preston at JPL [37]. 

Finally, the cycle-counting differential VLBI technique 
demonstrated in the ALSEP experiment has potential application 
for a number of future lunar and planetary missions [11] . 

For any mission in which transmitters on two or more space- 
craft are operating at nearly the same frequency, a DDR could 
be used for differential observations. If the spacecraft are 
equipped with transponders or stable oscillators, Doppler 
tracking may be done. In any case, differential interferometry 
may be used. The Viking project, in which two orbiters and 
two landers will be sent to Mars later this year, will provide 
several important applications for the techniques discussed in 

this thesis. 
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APPENDICES 

A. 1 Units and Pa r ameter Value s in the Theoretical Mode 1 
A . 1 . 1 Units 

The physical units of time, mass, and 
length used in the Planetary Ephemeris Program are de- 
fined as follows [ 6 ] : 

Time - The A.l second; defined as 9,192,631,770 
cycles at the terrestrial transition fre- 
quency of cesium at zero magnetic field 
[6 ]. 

Coordinate time (CT) , the independent 
variable of the equations of motion, is 
related theoretically to A.l time by 
CT = A.l + 32.15 sec + periodic general 
relativistic terms. 

The theoretical derivation of the rela- 
tivistic terms is discussed by Ash [ 6 ] . 

The definition presently implemented in 
PEP is 

CT - A.l + 32.15 sec + (1.658 x 10 ^sec) sin M 
+ (1.672 x 10 6 sec)sinD + i-~(r -r ) 

4 S G 
C 

where M is the earth's mean anomaly, D is 
the elongation of the moon from the sun, c 
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is the speed of light, r is the position of 

o 

the observing site with respect to the earth's 
center of mass, and r e is the velocity of the 
earth' s center of mass with respect to the sun. The 
coefficients of the annual and monthly terms are 
from Moyer [24 ] , and the expressions used for 
M and D are conventional [171* The vectors 
in the last (diurnal) term are evaluated using 
PEP's planetary ephemeris (NB0DY311) . 

Mass - The solar mass 

Length - The astronomical unit (AU) , defined as 

the radius of a heliocentric circular Kep- 
lerian orbit in which a body of negligible 
mass would have a period of 2ir/G days, 
where the gravitational constant, G, is as- 
sumed to be fixed at the Gaussian value, 

( .01720209895) 2 AU 3 /day 2 (1 day - 86400 A.l 
seconds). The accuracy of PEP's ephererides 
depends mainly on the measurements of the echo 
time delays of radar signals reflected from the 
inner planets. Simultaneously/ one also ob- 
tains the speed of light ( i . e . , of the radar 
signals) measured in AU/day; this result can 


be written conveniently as: 
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1 AU = 499.004780 ± 3 x 10~ 6 light seconds 
The conversion to metric units is accom- 
plished using the adopted value of the speed 
of light: 

c = 299792.458 km/sec 

A. 1.2 Parameters of the Lunar Orbit . The orbit 
of the moon was determined by fitting to laser ranging 
data [ 1, 33, 34] a numerical integration of the funda- 
mental equations of motion, as described by Ash [3 ] and 
Slade [36 ] . The parameter values determined from or used 
in this fit are as follows (not all digits given are 
necessarily significant) : 

Mass Ratio (sun)/ (earth-Hnoon) = 328900. 5011220281 

Mass ratio (earth/moon) = 81.3007 

Planetary masses and ephemerides - from NBODY311 

Osculating Keplerian elements at epoch JED 2440000.5 
(0 Hrs, CT on Julian Day Number 2440001) , with 
respect to 1950.0 equator and equinox: 

Semimajor axis, a = 0.002571513220080629 
= 0.05561587179127649 
= 28?396856 42497140 

= 3?31291959239597 
=226 ?271 133323205 


Eccentricity, e 

Inclination, i 

Right ascension 
of ascending 
node, ft 

Argument of 
perigee, oi 
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Initial mean 


anomaly, Mg 

= 154?885673245178 

Earth parameters: 


J 2 

= 1.082628 x 10" 3 

J 3 

=-2.538 x 10“ 6 

Radius used with 

harmonics, R 
e 

= 6378.166 km 

Tidal friction 
parameter, 
sin 26 

= 0.0755 (see Slade [36 ].) 

Moon parameters: 


J 2 

= 2.03822 x 10 -4 

C 22 

= 2.239570773993590 x 10" 5 

Radius used with 
harmonics, 

= 1738.09 km 

Inclination of 
lunar equator to 
ecliptic, I 

= 5552.71" 

Model of 


physical li- 
bration 

- Koziel, 1962 [20] 

1.3 Parameters of 

JPL Lunar-Libration Model 

(LLB-5 ) 


The libration 

model was generated by Wil- 

Sinclair of the Jet 

Propulsion Laboratory by 


fitting a numerical integration [40] to laser ranging data, 
The detailed characteristics of the fit are described in a 
memo [41 ] , available to the user . The nominal values of 
the parameters of the model are listed below for reference 


aiiitT-an- jit -ir 
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Moment of inertia ratios: 


3 = 631.26 x 10 
Y = 227.37 x 10' 


Inclination of lunar 
equator to ecliptic 


= 5552V 71 

“6 

Harmonics (units of 10 ) 

• 



J 2 

= 203.822 

C30 = - 10 - 44 

C 22 

= 22.3% 

C 31 = 28.6 

S 31 

= 8.8 

C32 = 4 ' 82 

S 32 

= 1.71 

C 33 = 2 ‘ 7 

S 33 

=-1.14 


Radius used with harmon- 
ics, R = 1738.09 km 

m 

Harmonics consistent with 
C/MR 2 = 0.394 

Mass of earth + moon = 328900.526 inverse 

solar mass units 

M as s ratio earth/moon = 81.3007 

Lunar ephemeris: LURE-2 


Initial conditions at epoch 
JED 2440400.5: 


p o = 
Ia o " 


% = 


1.245172332504406 x 10 _3 rad 

-3 

-0.070227621267468 x 10 rad 
0.1703911685703775 x 10 _3 rad 
-2.656265725254470 x 10~ 5 rad/day 
-1.948246963467498 x 10~ 5 rad/day 
-1.651226963560306 x 10 _5 rad/day 


r 


A. 2 List of Observations 



i 




The following list shows all of the differential 
VLB I observations of ALSEP transmitters which were made 
with MIT differential Doppler receivers at NASA Space- 
craft Tracking and Data Network (STDN) stations 
between Oct 1972 and June 1974. Date is the UTC date 
of the start of a continuous series of observations, 
that is, a series within which the differential Doppler 
counts are uninterrupted at both stations. The times 
are the UTC start and stop times; if the listed stop 
time of a series is numerically less than the start 
time , it is understood that the series stopped on the 
following UTC day. Stations are: GDS: the 85' "Apollo" 

station at Goldstone, California; TEX: the 30' station, 
now closed, at Corpus Christi, Texas; MIL: the original 
30' station on Merritt Island, Florida; ETC: the original 
30' station at Greenbelt, Maryland; MAD: the 85' station 
at Madrid, Spain; and ACN: the 30' station on Ascension 
Island in the South Atlantic Ocean. The first ALSEP 
listed is the "reference" ALSEP, that is, the one to 
whose carrier signal the S-band receiver was phase- 
locked. The second ALSEP is the one to whose carrier 
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signal the DDR was locked. The ALSEPs are designated 

by their Apollo mission numbers. A number^in the Obs. 

Lib. Tape column identifies the MIT observation library 

tape on which the data from a series is written. The 

« ' 

Series Name is formed from the day of the year and an iden- 
tifying letter (series in 1973 end in A through L, those 

in 1974 end in M through Z) . 

The Notes indicate whether or not the observation 
series was included in the data sets processed for this 
thesis. If a series was not included, the reason is indi- 
cated. The key to the notes is as follows: 

1 - The series was included in Data Set A. 

2 - The series was included in Data Sets B and C. 

3 — The series was included in Data Set C (as well 

as DS A) . 

4 - An error was made in preparing the data tapes to 
be used by PEP. The data are recoverable unless an error 
(not yet detected) of the type described by Note 5 also 
occurred. 

5 - Valid data were not obtained due to operator error 
or equipment malfunction. 

6 - All data taken at GDS between 22 Aug and 20 Sept 
1973 were corrupted by severe fluctuations (~ 1 part in 
10 10 ) in the primary frequency standard. 
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Times Stations ALSEPS Ob s.Li b^ Series. Note s 

uaxe Tape Name ;■ 

(mm/dd) (hhmm) ' y , 

1972 


10/28 

1973 

0634-1330 

MIL-GDS 

12-14 

10 

302A 

Omitted 

from 

these 

solutions 

03/7 

1350-2045 

MAD -ETC 

15-14 

11 

066G 

1 

03/8 

1352-2048 

1130-2048 

1352-2214 

1352-2048 

1130-2048 

1352-2214 

ACN-ETC 

acn-mad 

etc-mad 

ACN-ETC 

acn-mad 

etc-mad 

15-14 

15-14 

15-14 

15-12 

15-12 

15-12 

11 

11 

11 

12 

12 

12 

067G 
067H 
0671 
067 J 
067K 
067L 

1 

1 

4 

1 

1 

1 

03/15 

2038-0335 

0146-0330 

0146-0330 

2045-0335 

2045-0335 

2045-0335 

MAD -MIL 
MAD-ACN 
MIL-ACN 
MAD -MIL 
MAD-ACN 
MIL-ACN 

15-14 
15-14 
15-14 
IS- 12 
15-12 
15-12 

11 

11 

11 

12 

12 

12 

074G 

075H 

0751 

075J 

075K 

075L 

4 

5 
5 
5 
1 
5 

03/22 

2348-0404 

ACN-MAD 

15-14 

11 

08 IB 

1 

03/24 

0215-0545 

0215-0447 

ACN-MAD 

ACN-MAD 

15-14 

15-17 

11 

13 

083B 

083A 

5 

5 

03/28 

1115-1915 

gds-tex 

15-14 

11 

087B 

1 

03/30 

1135-1415 

1135-1415 

ACN-MAD 

acn-mad 

15-14 

15-17 

11 

13 

089D 

089B 

5 

5 

04/3 

1336-1639 

tex-mad 

15-14 

11 

093G 

4 

04/4 

1454-1627 

mad-tex 

15-14 

11 

094G 

4 

04/9 

1815-2330 

2351-0400 

ACN-TEX 

etC-tex 

15-17 

15-17 

13 

13 

099A 

099B 

5 

5 

04/12 

2115-0059 

2110-0059 

MAD-ACN 

MAD-TEX 

15-14 

15-14 

11 

11 

102A 

102B 

4 

4 

04/13 

1925-0245 

1635-0245 

acn-mad 

ACN-MAD 

15-14 

15-16 

49 

48 

10 3 A 
103B 

1 

5 
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Date 
(mm/ dd ) 

Times 

(hhmm) 

Stations 

ALSEPS 

Obs. Lib. 
Tape 

Series 

Name 

04/17 

2040-0300 

ACN-MAD 

15-14 

11 

107A 

0228-0520 

MAD -TEX 

15-14 

11 

107B 


0435-1215 

TEX-GDS 

15-14 

11 

107C 

04/18 

2100-0300 

ACN-MAD 

15-14 

11 

108D 

2100-0300 

ACN-MAD 

15-17 

13 

108A 


0230-0545 

MAD-TEX 

15-17 

13 

108B 


0435-1300 

TEX-GDS 

15-17 

13 

108C 

07/2 

1450-2035 

ACN-TEX 

15-14 

21 

183A 

1600-2030 

TEX-GDS 

15-14 

21 

183B 


2035-2335 

TEX-GDS 

15-14 

21 

183D 


2335-0305 

TEX-GDS 

15-14 

21 

183E 


2115-2330 

GDS— MIL 

15-17 

22 

183C 

07/3 

1815-0200 

MIL-GDS 

15-14 

21 

184A 

2015-0200 

MIL-TEX 

15-14 

21 

184B 


1530-2000 

TEX-MAD 

15-14 

21 

184C 


2005-0345 

TEX-GDS 

15-14 

21 

184D 


1815-0200 

MIL-GDS 

15-17 

22 

184E 

07/5 

1230-1600 

MAD-ACN 

15-17 

22 

186A 

1955-0315 

MIL-GDS 

15-12 

24 

18 6B 


0000-0500 

GDS-TEX 

15-16 

23 

186C 


1230-1700 

MAD-AON 

15-16 

25 

186D 


1830-2200 

MAD-ACN 

15-16 

25 

186H 


1645-2200 

mad-mil 

15-16 

25 

186E 


1955-0315 

MIL-GDS 

15-16 

25 

186F 


1800-2315 

MIL-TEX 

15-16 

25 

186G 

07/6 

1410-2140 

ACN-MAD 

14-15 

21 

187A 

1715-0015 

ACN-MIL 

14-15 

21 

187B 


2050-0245 

MIL-GDS 

14-15 

21 

187C 


1410-2140 

ACN-MAD 

14-12 

23 

187D 


1715-0015 

ACN-MIL 

14-12 

23 

187E 


2050-0245 

MIL-GDS 

14-12 

23 

187F 

07/9 

2145-0245 

ACN-TEX 

15-16 

25 

190A 

2225-0315 

MIL-TEX 

15-16 

25 

19 OB 


0325-0730 

TEX-GDS 

15-16 

25 

190C 


0245-0545 

MIL-GDS 

15-17 

22 

190D 


2225-0235 

ACN-MIL 

15-17 

22 

190E 


Notes 


5 

5 

5 


5 

5 

1 

1 

1 

1 

1 

1 

5 

1 

1 

5 

1 

4 

4 

4 
1 
1 

5 
1 
1 
1 


<sf ia h h 
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Date. 

(mm/dd) 

07/10 

07/11 

07/12 

07/13 

07/17 

07/18 

07/19 


Times 

Stations 

AI.SEPS 

(hhmm) 

2100-0115 

MAD-MIL 

15-14 

0015-0330 

MIL-GDS 

15-14 

2100-0115 

MAD-MIL 

15-16 

0015 **06 00 

MIL-GDS 

15-16 

0015-0600 

TEX-GDS 

15-16 

2200-0145 

MAD-MIL 

15-17 , 

0200-0600 

MIL-TEX 

15-17 

0115-0600 

MIL-GDS 

15-17 

0115-0600 

MIL-GDS 

15-12 

1845-0050 

ACN-MAD 

15-16 

2245-0515 

ACN-MIL 

15-16 

0200-0630 

MIL-GDS 

15-16 

0000-0430 

MIL-TEX 

15-16 

0115-0515 

ACN-MIL 

15-12 

1845-0050 

ACN-MAD 

15-12 

0150-0630 

MIL-GDS 

15-12 

1946-0330 

ACN-MAD 

15-17 

2355-0300 

ACN-MIL 

15-17 

0240-0553 

ACN-GDS 

15-17 

1946-0330 

ACN-MAD 

15-12 

2355-0300 

ACN-MIL 

15-12 

0240-0600 

ACN-GDS 

15-12 

0220-0630 

MAD-ACN 

15-17 

0230-0730 

ACN-TEX 

15-17 

0415-0730 

ACN-GDS 

15-17 

0415-1030 

TEX-GDS 

15-17 

0220-0630 

MAD-ACN 

15-16 

0415-0730 

ACN-GDS 

15-16 

0315-0735 

TEX-MAD 

15-12 

0452-1415 

TEX-MAD 

15-12 

0500-0735 

GDS-MAD 

15-17 

0345-0810 

MAD-TEX 

15-12 

0530-1205 

TEX-GDS 

15-12 

0530-0810 

MAD-GDS 

15-16 


Obs. Lib. 

Series 

Note: 

Tape 

Name 


11 

191A 

5 

11 

191B 

5 

25 

191C 

4 

25 

191D 

1 

25 

191E 

4 

22 

192A 

4 

22 

192B 

4 

22 

192C 

4 

24 

192D 

4 

25 

193A 

1 

25 

193 B 

1 

25 

193C 

5 

25 

19 3D 

4 

24 

193E 

4 

24 

193F 

4 

24 

193G 


22 

194A 

1 

22 

194B 

4 

22 

194C 

4 

24 

194D 

4 

24 

194E 

4 

24 

194F 

4 

22 

198A 

1 

22 

198B 

1 

22 

198C 

1 

22 

198D 

4 

25 

198E 

1 

25 

198F 

5 

24 

199A 

4 

24. 

199B 

4 

22 

199C 

4 

24 

200A 

4 

24 

200B 

4 

25 

200C 

4 
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Date Times stations 

(mm/dd) (hhmm) 


07/20 

0000-0700 

MAD-ACN 

0620-0930 

ACN-GDS 


0620-0930 

ACN-GDS 


0000-0700 

MAD-ACN 

07/23 

0615-1150 

TEX-ACN 

1425-1850 

GDS-TEX 


0930-1225 

ACN-GDS 

07/24 

0230-0730 

ACN-MAD 

0820-1400 

MAD-GDS 


0820-1745 

TEX-GDS 


0230-0730 

ACN-MAD 


0820-1400 

MAD-GDS 

07/25 

0745-0912 

MAD-TEX 

0850-1700 

GDS-TEX 


1000-1515 

MAD-GDS 

08/6 

2315-0415 

GDS-MIL 

2318-0415 

TEX-GDS 


2315-0415 

GDS-MIL 

08/8 

0015-0500 

MIL-GDS 

0015-0500 

TEX-GDS 


0015-0500 

MIL-GDS 


0000-0330 

MIL-GDS 


2213-0330 

MIL-TEX 


0000-0330 

TEX-GDS 

08/9 

1800-2055 

MAD-ACN 

0055-0630 

MIL-GDS 


2330-0500 

MIL-TEX 


0240-0630 

MIL-GDS 


1800-2055 

MAD-ACN 

08/14 

0105-0600 

MIL-ACN 

0105-0600 

MIL-ACN 

08/17 

0240-0700 

MIL-ACN 

0240-0700 

MIL-ACN 

08/21 

0600-1305 

MAD -MIL 

1435-1730 

MIL-GDS 


0600-1100 

mad-mil 


1405-1800 

GDS-TEX 


1325-1730 

MIL-TEX 


ALSEPS 

Obs . Lib . 

Series 

Notes 


Tape 

Name 


15-14 

49 

201A 

4 

15-14 

49 

201B 

4 

15-16 

48 

201C 

4 

15-16 

48 

201D 

4 

12-17 

26 

204A 

4 

12-17 

26 

204B 

2 

12-14 

21 

204C 

2 

12-17 

26 

205A 

2 

12-17 

26 

205B 

4 

12-17 

26 

205C 

2 

12-14 

23 

205D 

4 

12-14 

23 

205E 

4 

12-17 

26 

206A 

5 

12-17 

26 

206C 

2 

12-14 

23 

206D 

2 

12-14 

23 

2 18 A 

4 

12-14 

23 

218B 

4 

12-17 

26 

218C 

5 

15-17 

22 

220A 

1 

15-17 

22 

2 2 OB 

1 

15-16 

25 

220C 

4 

15-17 

22 

220D . 

4 

15-16 

25 

220E 

5 

15-16 

25 

220F 

5 

12-17 

26 

221A 

4 

12-17 

26 

221B 

2 

12-15 

24 

221C 

1 

12-15 

24 

221F 

4 

12-15 

24 

2 2 IE 

4 

12-15 

24 

226A 

4 

12-14 

23 

226B 

5 

12-15 

24 

229A 

4 

12-14 

23 

229B 

4 

12-14 

23 

233A 

4 

12-15 

24 

233B 

6 

12-17 

26 

233C 

4 

12-17 

26 

233D 

6 

12-17 

26 

233E 

4 
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Date 

Times 

Stations 

ALSEPS 

Obs. Lib. 

Series 

Notes 

(nun/dd) 

(hhmm) 



Tape 

Name 


08/22 

0603-1400 

MAD-MIL 

16-17 

47 

234A 

4 

1300-1830 

MIL-GDS 

16-17 

47 

234B 

6 


0603-1400 

MAD-MIL 

16-15 

25 

234C 

1 


1300-1830 

MIL-GDS 

16-15 

25 

234D 

6 


1316-2000 

GDS-TEX 

16-15 

25 

234E 

6 

08/23 

1440-1620 

TEX-GDS 

12-17 

26 

235A 

6 

1700-2115 

TEX-GDS 

12-17 

26 

235B 

6 

08/24 

0932-1440 

ACN-MIL 

12-15 

24 

236A 

5 

1300-1945 

MIL-GDS 

12-15 

24 

236B 

6 


0932-1440 

ACN-MIL 

12-14 

23 

236C 

2 


1300-1945 

MIL-GDS 

12-14 

23 

236D 

6 

08/27 

1502-2200 

MIL-GDS 

12-14 

23 

239A 

6 

1502-2200 

MIL-GDS 

12-17 

26 

239B 

6 


1502-2028 

GDS-TEX 

12-17 

26 

239C 

6 

08/28 

1508-2300 

MIL-GDS 

12-14 

23 

240A 

6' 

1508-2300 

MIL-GDS 

12-15 

24 

240B 

6 


1508-2130 

TEX-GDS 

12-15 

24 

240C 

6 

08/29 

1600-2230 

MIL-GDS 

12-14 

23 

241A 

6' 

1600-2230 

MIL-GDS 

12-15 

24 

241B 

6 


1600-2100 

TEX-GDS 

12-15 

24 

241C 

6 

08/30 

1633-2130 

MIL-GDS 

12-17 

26 

242A 

6 

1636-2130 

MIL-GDS 

12-16 

34 

242B 

6 


1508-2000 

MIL-TEX 

12-16 

34 

242C 

2 

09/3 

1830-2400 

ACN-MIL 

12-15 

24 

246A 

1 

1830-2400 

ACN-MIL 

12-14 

23 

246B 

2 

09/10 

0005-0400 

MIL-MAD 

15-16 

48 

253A 

1,3 

2352-0400 

MIL-MAD 

15-14 

49 

253B 

5 

09/11 

2345-0450 

MAD-MIL 

15-16 

48 

254A 

1, 3 

0600-1015 

MIL-GDS 

15-16 

48 

254B 

6 


0605-1200 

TEX-GDS 

15-16 

48 

254C 

6 


2345-0445 

MAD-MIL 

15-17 

46 

254D 

1 


0600-1015 

MIL-GDS 

15-17 

46 

254E 

6 

09/13 

0015-0530 

MAD-MIL 

15-17 

46 

256A 

1 

0500-1030 

MIL-TEX 

15-17 

46 

256B 

1 


0630-1255 

GDS-TEX 

15-17 

46 

256C 

6 


0015-0530 

MAD-MIL 

15-14 

49 

256D 

1 


0620-1030 

MIL-GDS 

15-14 

49 

256E 

6 
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Date 

Times 

Stations 

ALSEPS 

Obs. Lib. 
Tape 

Series 

Name 

Notes 

(mm/dd) 

(hhmm) 





09/14 

0125-0550 

0034-0630 

0605-1200 

0125-0550 

0100-0630 

0605-1200 

MAD- AC N 
MIL-MAD 
MIL-GDS 
MAD-ACN 
MIL-MAD 
MIL-GDS 

12-15 

12-15 

12-15 

12-14 

12-14 

12-14 

32 

32 

32 

33 
33 
33 

257A 

257B 

257C 

257D 

257E 

257F 

1 

1 

6 

2 

2 

6 

09/15 

0115-0700 

0115-0700 

ACN-MAD 

acn-mad 

12-16 

12-17 

34 

31 

258A 

258B 

4 

5 

09/17 

0640-1200 

0858-1330 

0858-1200 

MIL-TEX 

gds-tex 

MIL-GDS 

12-15 

12-15 

12-14 

32 

32 

33 

260A 

260B 

260C 

5 

6 
6 

09/18 

0510-1200 

MIL-TEX 

12-17 

31 

261A 

2 

09/19 

0540-0930 

0935-1500 

0630-1500 

MIL-TEX 

MIL-GDS 

MIL-GDS 

12-17 

12-17 

12-14 

31 

31 

33 

262A 

262B 

262C 

4 

6 

6 

09/20 

1126-1340 

1013-1230 

1235-1430 

MIL-GDS 

MIL-GDS 

MIL-GDS 

12-15 

12-16 

12-16 

32 

34 

34 

263A 

263B 

263C 

6 

6 

6 

09/24 

0930-1530 

0930-1530 

acn-mil 

ACN-MIL 

12-15 

12-17 

32 

31 

267A 

267B 

1 . 

4 

09/26 

1230-1645 

1352-2200 

1230-1645 

MAD-MIL 

MIL-TEX 

MAD-MIL 

12-17 

12-17 

12-15 

31 

31 

32 

269A 

269B 

269C 

4 

2 

1 

09/27 

0830-1730 

1350-1900 

0830-1730 

1350-1900 

acn-mad 

acn-mil 

acn-mad 

acn-mil 

12-17 

12-17 

12-16 

12-16 

31 

31 

34 

34 

270A 

270B 

270C 

270D 

2 

2 

2 

2 

09/28 

0916-1800 

1402-2000 

0916-1800 

1402-2000 

acn-mad 

acn-mil 

acn-mad 

acn-mil 

12-16 

12-16 

12-14 

12-14 

34 

34 

33 

33 

271A 

271B 

271C 

271D 

2 

4 

2 

2 

10/1 

1720-2245 

2235-0130 

1720-2245 

acn-mil 

MIL-TEX 

ACN-MIL 

12-16 

12-16 

12-15 

34 

34 

32 

274A 

274B 

274C 

4 

2 

1 

10/2 

1720-2100 

2100-0230 

1720-2100 

MIL-MAD 

MIL-TEX 

MIL-MAD 

12-15 

12-15 

12-16 

32 

32 

34 

275A 

275B 

275C 

1 

1 

2 
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Date 

Times 

Stations 

ALSEPS 

Obs. Lib. 

Series 

Notes 





Tape 

Name 


(mm/dd) 

(hhmm) 






10/3 

1820-2200 

MIL-MAD 

12-15 

32 

276A 

1 

T 


2010-0300 

MIL-TEX 

12-15 

32 

276B 

X 


1835-2200 

MIL-MAD 

12-14 

33 

276C 

4 

10/4 

1930-2300 

ACN-MAD 

12-16 

34 

277A 

4 

A 


1930-2300 

ACN-MIL 

12-16 

34 

277B 



0025-0400 

MIL-TEX 

12-16 

34 

277C 

D 


1930-2300 

ACN-MAD 

12-14 

33 

277D 

Z 


1930-2300 

ACN-MIL 

12-14 

33 

2 77E 

z 

10/5 

1730-2200 

ACN-MAD 

12-17 

31 

278A 

2 


1923-2400 

MAD-MIL 

12-17 

31 

278B 

5 


1730-2200 

ACN-MAD 

12-14 

33 

278C 

2 


1923-2400 

MAD-MIL 

12-14 

33 

278D 

4 

10/16 

0300-0700 

ACN-MIL 

12-15 

32 

289A 

1 


0300-0700 

ACN-MIL 

12-14 

33 

289B 

2 

10/22 

0950-1430 

ACN-MIL 

12-17 

31 

295A 

4 


0950-1430 

ACN-MIL 

12-16 

34 

295B 

2 

10/25 

1140-1600 

MAD-MIL 

12-15 

32 

298A 

1 


1140-1600 

MAD-MIL 

12-14 

33 

298B 

2 

11/2 

1500-2200 

MAD-ACN 

12-17 

31 

306A 

4 


1905-2230 

MAD-MIL 

12-17 

31 

306B 

4 


1500-2200 

MAD-ACN 

12-14 

33 

306C 

2 


1905-2230 

MAD-MIL 

12-14 

33 

306D 

2 

11/5 

2320-0600 

MIL-GDS 

16-17 

47 

309A 

2 


2325-0600 

MIL-GDS 

16-15 

48 

309B 

1 1 3 

11/7 

0000-0700 

mil-gds 

16-17 

47 

311A 

2 


0000-0700 

MIL-GDS 

16-15 

48 

311B 

1 1 3 


2130-0300 

MAD-MIL 

12-17 

31 

311C 

4 


0030-0550 

MIL-GDS 

12-17 

31 

3 11D 

4 


2130-0300 

MAD-MIL 

12-14 

33 

311E 

4 


0030-0800 

MIL-GDS 

12-14 

33 

31 IF 

2 

11/8 

1832-0415 

MAD-ACN 

12-15 

32 

312A 

4 


2150-0500 

MAD-MIL 

12-15 

32 

312B 

4 


0030-0700 

MIL-GDS 

12-15 

32 

312C 

4 


1832-0415 

MAD-ACN 

12-14 

33 

312D 

4 


2150-0500 

MAD-MIL 

12-14 

33 

312E 

2 


0030-0700 

MIL-GDS 

12-14 

33 

312F 

2 
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Date 

(mm/dd) 

Times 

(hhmm) 

Stations 

ALSEPS 

Obs . Lib . 
Tape 

Series 

Name 

Notes 

11/9 

1934-0500 

2202-0400 

0045-0600 

1945-0500 

2202-0400 

0045-0600 

MAD-ACN 
MAD -MIL 
MAD-GDS 
MAD-ACN 
MAD -MIL 
MAD-GDS 

12-16 

12-16 

12-16 

12-15 

12-15 

12-15 

34 

34 

34 

32 

32 

32 

313A 

313B 

313C 

313D 

313E 

313F 

4 

2 

2 

4 

1 

1 

11/27 

0500-1830 

1500-1830 

MAD-MIL 
MAD -MIL 

12-15 

12-14 

44 

42 

331A 

331B 

1 

2 

11/28 

1545-1945 

1545-1940 

MAD-MIL 

MAD-MIL 

12-17 

12-16 

41 

43 

332A 

332B 

2 

2 

11/29 

1700-2000 

1700-2200 

2155-0115 

1756-2000 

1756-2115 

2155-0115 

MAD-ACN 

MIL-ACN 

MIL-GDS 

MAD-ACN 

MIL-ACN 

mil-gds 

12-17 

12-17 

12-17 

12-16 

12-16 

12-16 

41 

41 

41 

43 

43 

43 

333A 

333B 

333C 

333D 

333E 

333F 

5 

5 

2 

5 

5 

5 

11/30 

1750-2130 

2145-0230 

1745-2130 

1715-2130 

2145-0230 

MAD-MIL 

MIL-GDS 

ACN-MAD 

MAD-MIL 

MIL-GDS 

12-15 

12-15 

12-14 

12-14 

12-14 

44 

44 

42 

42 

42 

334A 

334B 

334C 

334D 

334E 

1 

1 

4 

5 
5 

12/1 

1645-2100 

1645-2100 

ACN-MAD 

ACN-MAD 

12-15 

12-14 

44 

42 

335A 

335B 

1 

2 

12/3 

2305-0400 

2305-0400 

MIL-GDS 

MIL-GDS 

15-17 

15-16 

46 

45 

337A 

337B 

4 

4 

12/4 

1915-0030 

2315-0600 

1915-0030 

2309-0600 

MIL-MAD 

MIL-GDS 

MIL-MAD 

MIL-GDS 

15-17 

15-17 

15-16 

15-16 

46 

46 

45 

45 

338A 

338B 

338C 

338D 

4 

4 

4 

4 

12/6 

0218-0650 

0218-0650 

MIL-GDS 

MIL-GDS 

15-16 

15-16 

45 

45 

340A 

340B 

4 

If 3 
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Date 

Times 

Stations 

(mm/dd) 

(hhitun) 


12/7 

0230-0645 

0230-0645 

MIL-GDS 

MIL-GDS 

12/20 

1426-1830 

1426-1830 

MIL-GDS 

MIL-GDS 

12/21 

1330-1915 

1321-1915 

MIL-GDS 

MIL-GDS 

1974 



01/- 

0345-1515 

0345-1515 

0530-1200 

0530-1200 

MAD-MIL 

MAD-MIL 

MIL-GDS 

MIL-GDS 

01/10 

0215-0600 

0455-1100 

0215-0600 

0500-1100 

MAD-MIL 

MIL-GDS 

MAD-MIL 

MIL-GDS 

01/11 

0200-0600 

0605-1100 

0200-0600 

0605-1100 

MAD-MIL 

MIL-GDS 

MAD-MIL 

MIL-GDS 

01/28 

2215-0230 

2215-0230 

MIL-GDS 

MIL-GDS 

01/30 

1800-2400 

2112-0400 

1800-2400 

2112-0400 

MAD-MIL 

MIL-GDS 

MAD-MIL 

MIL-GDS 

02/1 

1850-2400 

2100-0145 

2100-0145 

MAD-MIL 

MIL-GDS 

MIL-GDS 

02/2 

1746-2200 

1900-2200 

MAD -AC N 
MAD-ACN 

02/15 

0515-1000 

0447-1000 

ACN-MAD 

ACN-MAD 


ALSEPS 

Obs. Lib. 

Series 

Notes 


Tape 

Name 


12-17 

41 

340C 

4 

12-14 

42! 

340D 

4 

12-17 

41 

354A 

5 

12-14 

42 

354B 

2 

12-15 

44 

355A 

1 

12-16 

43 

355B 

2 


12-17 

41 

009M 

12-17 

41 

009N 

12-17 

41 

009P 

12-17 

41 

009Q 

15-16 

45 

010M 

15-16 

45 

010N 

15-12 

44 

010P 

15-12 

44 

010Q 

14-12 

42 

OHM 

14-12 

42 

011N 

14-12 

42 

011P 

14-12 

42 

011Q 

12-16 

43 

028M 

12-14 

42 

028N 

12-17 

41 

030M 

12-17 

41 

030N 

12-14 

42 

030P 

12-14 

42 

030Q 

16-17 

47 

032M 

16-17 

47 

032N 

16-15 

45 

032P 

16-17 

47 

033M 

16-15 

45 

033N 

12-17 

41 

046M 

12-16 

43 

046N 


02/16 


0430-1000 ACN-MAD 15-17 46 

0425-1000 ACN-MAD 15-16 45 


047M 

047N 


5 

4 


*,*,*... NJ NJ ^ U1 HH' - 
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Date 

Times 

(mm/dd) 

(hhmm) 

02/20 

1405-1710 

02/21 

1015-1700 

1400-1800 

1500-2200 

1135-1700 

1352-1800 

1454-2200 

02/22 

1200-1745 

1400-1800 

1509-2200 

1200-1745 

1400-1800 

1500-2200 

02/23 

1200-1830 

1200-1830 

02/25 

1348-2030 

1700-0100 

1348-2030 

1700-0100 

02/26 

2317-0200 

2317-0200 

02/27 

1745-0300 

1745-0300 

02/28 

1615-2007 

2010-0300 

1615-2100 

1824-0300 

03/4 

2235-0300 

2235-0300 

03/15 

0240-0745 

0240-0910 


Stations 

ALSEPS 

MIL-GDS 

12-15 

MAD-ACN 

12-17 

ACN-MIL 

12-17 

MIL-GDS 

12-17 

MAD-ACN 

12-16 

ACN-MIL 

12-16 

MIL-GDS 

12-16 

ACN-MAD 

12-15 

ACN-MIL 

12-15 

MIL-GDS 

12-15 

ACN-MAD 

12-14 

ACN-MIL 

12-14 

MIL-GDS 

12-14 

ACN-MAD 

12-17 

ACN-MAD 

12-14 

ACN-MIL 

12-17 

MIL-GDS 

12-17 

ACN-MIL 

12-16 

MIL-GDS 

12-16 

MIL-GDS 

12-15 

MIL-GDS 

12-16 

MIL-GDS 

12-15 

MIL-GDS 

12-14 

ACN-MIL 

12-17 

MIL-GDS 

12-17 

ACN-MIL 

12-14 

MIL-GDS 

12-14 

MIL-GDS 

15-17 

MIL-GDS 

15-16 

ACN-MAD 

12-17 

ACN-MAD 

12-14 


Obs . Lib . 

Series 

Tape 

Name 

44 

051N 

41 

052M 

41 

052N 

41 

052P 

43 

052Q 

43 

052R 

43 

052S 

44 

053M 

44 

053N 

44 

053P 

42 

053Q 

42 

053R 

42 

053S 

41 

054M 

42 

054N 

41 

056M 

41 

056N 

43 

056P 

43 

056Q 

44 

057M 

43 

057N 

44 

058M 

42 

058N 

41 

059M 

41 

059N 

42 

059P 

42 

059Q 

46 

063M 

45 

063N 

50 

074M 

51 

074N 


Notes 


5 


5 

1,3 

5 

2 


mcMCM^oicsj i-t r- i n* cn cm <n in in cm cm cm cm ^ *H cm in cm ^ cm 
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Date 

(rnm/dd) 

Times 

(hhmm) 

Stations 

03/22 

1300-1630 

1300-1630 

MAD -MIL 
MAD-MIL 

03/26 

0955-1425 

1315-1715 

0955-1425 

1315-1715 

ACN-MAD 

MAD-MIL 

ACN-MAD 

MAD-MIL 

03/27 

1720-2200 

1110-1630 

1400-2200 

MAD-MIL 

MAD-ACN 

MAD-MIL 

04/1 

2136-0200 

2040-0200 

2136-0200 

2040-0200 

ACN-MAD 

ACN-MIL 

ACN-MAD 

ACN-MIL 

04/2 

2010-0200 

2010-0200 

2010-0200 

2010-0200 

MAD-ACN 

MAD-MIL 

MAD-ACN 

MAD-MIL 

04/9 

0000-0630 

0000-0630 

2300-0415 

2315-0415 

ACN-MAD 

ACN-MAD 

ACN-MAD 

ACN-MAD 

04/16 

0400-1230 

0400-1230 

MAD-ACN 

MAD-ACN 

04/22 

1320-1845 

1250-1845 

MAD-MIL 

MAD-MIL 

04/23 

1445-2330 

1445-2330 

GDS-MIL 

GDS-MIL 

04/24 

1500-1930 

1455-1904 

ACN-MAD 

ACN-MAD 

04/29 

1845-2312 

1845-2400 

1845-2312 

1845-2400 

ACN-MAD 

MAD-MIL 

ACN-MAD 

MAD-MIL 


ALSEPS Obs. Lib. 

Series Notes 

Tape 

Name 


12-15 

52 

081M 

1 

12-16 

53 

08 IP 

2 

16-17 

54 

085M 

2 

16-17 

54 

085N 

2 

16-15 

55 

085P 

1/3 

16-15 

55 

085Q 

1/3 

15-14 

56 

086M 

1 

15-16 

55 

086N 

1/3 

15-16 

55 

086P 

1/3 

16-17 

54 

091M 

5 

16-17 

54 

091N 

5 

16-15 

55 

091P 

4 

16-15 

55 

091Q 

1/3 

16-17 

54 

092M 

5 

16-17 

54 

092N 

2 

16-15 

55 

092P 

1/3 

16-15 

55 

092Q 

1/3 

14-12 

51 

099M 

4 

14-15 

56 

099N 

4 

12-14 

51 

099P 

2 

12-17 

50 

099Q 

5 

17-16 

54 

106M 

2 

17-15 

57 

106N 

4 

17-15 

57 

112M 

5 

17-16 

54 

112N 

5 

15-12 

52 

113M 

1 

15-14 

56 

113N 

1 

17-15 

57 

114M 

1 

17-16 

54 

114N 

2 

16-12 

53 

119M 

5 

16-12 

53 

119N 

5 

16-17 

54 

119P 

5 

16-17 

54 

119Q 

5 
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Date 

Times 

Stations 

ALSEPS 

Obs. Lib. 

Series 

Notes 





Tape 

Name 


(mm/dd) 

(hhmm) 






.04/30 

1715-2300 

ACN-MAD 

16-12 

53 

120M 

4 

1715-2300 

ACN-MAD 

16-17 

54 

120N 

4 

05/1 

1955-2400 

MAD-MIL 

16-15 

55 

121M 

5 


1955-2400 

MAD -MIL 

16-17 

54 

12 IN 

5 

05/11 

0128-0800 

ACN-MAD 

14-12 

51 

13 1M 

5 


0128-0800 

ACN-MAD 

14-15 

56 

131N 

5 

05/12 

0145-0800 

ACN-MAD 

12-15 

52 

132M 

5 


0150-0800 

ACN-MAD 

12-17 

50 

132N 

5 

05/13 

0126-0700 

ACN-MAD 

12-14 

51 

133M 

5 


0126-0700 

ACN-MAD 

12-17 

50 

133N 

5 

05/15 

0228-0850 

ACN-MAD 

16-15 

55 

135M 

5 


0228-0850 

ACN-MAD 

16-17 

54 

135N 

5 

05/17 

1300-1845 

GDS-MIL 

12-15 

52 

137M 

5 


1627-1845 

GDS-MIL 

12-17 

50 

137N 

5 

05/20 

1155-1600 

MAD-ACN 

16-17 

54 

140M 

5 

1830-2200 

MIL-GDS 

16-17 

54 

140N 

5 


1155-1600 

MAD-ACN 

16-12 

53 

140P 

4 


1830-2200 

MIL-GDS 

16-12 

53 

140Q 

2 

05/21 

1908-2315 

MIL-GDS 

15-16 

55 

141M 

5 

1908-2315 

MIL-GDS 

15-14 

56 

141N 

1 

05/22 

1115-1800 

MAD-ACN 

15-16 

55 

142M 

4 

1115-1800 

MAD-ACN 

15-14 

56 

142N 

1 

05/27 

1240-2000 

MAD-ACN 

12-14 

51 

147M 

5 

1245-2000 

MAD-ACN 

12-15 

52 

147N 

1 

05/28 

1830-2200 

MAD-MIL 

12-14 

51 

148M 

5 

1830-2200 

MAD-MIL 

12-16 

53 

148N 

5 

05/30 

2100-0100 

MAD-MIL 

16-15 

55 

150M 

1,3 

2210-0210 

MIL-ACN 

16-15 

55 

150N 

1,3 


2100-0100 

MAD-MIL 

16-17 

54 

150P 

5 


2210-0210 

MIL-ACN 

16-17 

54 

150Q 

2 


Date 

Times 

Stations 

ALSEPS 

Obs. Lib. 

Series 





Tape 

Name 

(mm/dd) 

(hhmm) 





06/5 

0110-0400 

ACN-MAD 

16-15 

55 

156M 


0110-0400 

ACN-MAD 

15-15 

55 

156N 

06/7 

0005-0600 

ACN-MAD 

16-15 

55 

158M 


0005-0600 

ACN-MAD 

16-12 

53 

158N 

06/11 

0215-0935 

MAD-ACN 

12-15 

52 

162M 


0215-0935 

MAD -AC N 

12-14 

51 

162N 

06/13 

0215-1115 

ACN-MAD 

12-17 

50 

164M 


0215-1115 

ACN-MAD 

12-16 

53 

164N 

06/17 

1300-2000 

MIL-GDS 

16-17 

54 

168M 


1300-2000 

MIL-GDS 

16-12 

53 

168N 

06/18 

0725-1430 

MAD-ACN 

16-17 

54 

169M 


0725-1430 

MAD-ACN 

16-12 

53 

169N 

06/19 

0700-1335 

MAD-ACN 

15-14 

56 

17 0M 


0700-1330 

MAD-ACN 

15-17 

57 

170N 

06/20 

1830-2230 

MIL-GDS 

15-14 

56 

171M 


1830-2230 

MIL-GDS 

15-17 

57 

171N 


Notes 


5 

4 

4 

2 


5 

2 

5 

2 

4 

4 

1 

4 

5 
1 


.H (N 


A. 3 


Derivation of Terms in the Differential N-Count 
Observable Resultin g from the Moon's Horizontal 
Parallax 

This appendix presents the derivation of Equation (2.2.8) 
which includes the effect of horizontal parallax in the dif- 
ferential N-count observable for the case where both the 
tracking stations and the ALSEPs lie in the plane of the 
earth's equator. Figure A. 3.1 shows the geometry for a 
single station observing a single transmitter on the moon. 



Figure A. 3.1 
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In the figure, R is the earth-moon center-to-center 
distance, r Q is the earth's radius, r m is the moon's 
radius, is the selenocentric longitude of the trans- 
mitter, and hj is the hour angle of the moon at the track- 
ing station. The differential fringe phase, or differential 
N-count, may be written 


At}) (t) = — { [^11 (t) ~ (t:) ] ™ [D 9 ^ (t) - 022 (b)]} (A. 3.1) 

where (t) is the distance from transmitter i to 
station j. Here, and in the remainder of this Appendix, 
the constant term A4> , which accounts for the initial 
counter readings is omitted for simplicity. 

D 11 = ( R " r g cos h 1 - r^cos Jl^) 2 + ( r m sin + r e sin h i^ 2 

= R 2 + r 2 cos 2 h 1 + r 2 cos 2 £ 1 - 2Rr e cos - 2Rr m cos 

2 2 2.2 

+ 2r r cos h.cos + r sin + r sin h. 
em 1 l m 1 e i. 

+ 2r r sin S,, sin h-, (A. 3.2) 

me 1 l 


By combining terms and using the identities 

sin 2 x + cos 2 x = 1 and cos (x - y) = cos x cos y + sin x sin y. 
Equation (A. 3.2) can be written 

0 2 1 = r2[1 + (~) 2 + (“) 2 - 2(jp-)cos h 1 -2 (”) cos SL 1 


r r 

+ 2 ( %~ )cos(Jl, - h. ) ] 
IT 


(A. 3. 3) 
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Since <=£) « .017 and <^> * .004, the expression in 
R *'■ 


R ' R 

brackets in Equation (A. 3. 3) has the form (1 + x> with 

1/2 _ 


|x| << 1, and we can use the expansion (1 + x) 


! + | x - i-x 2 + jlx 3 ... for the square root of D^. Ex- 


pending, and retaining all terms through order (r/R) , where 


r is r e or r^, we obtain 


Dll = ru * |(^> 2 - |<^> 2 ' < ir> c ° s h i - ( ^ )cos 


1 , r nu 2 


+ (Is^DcostVh,) - |(^) 2 cos 2 h 1 - 

R 2 r 

- ( e Jl)cos h, cos 2^ + (-—--} cos h 1 cos^-h^) 
R 2 R 


+ (I®™) cos l 1 cos( 2 - 1 -h 1 ) + cos h i 

R J 


2 2 

+ |<^) 3 cos 4 + |(^)cos h, + i(^?)cos h 


R’ 


y* IT 

1 r e 3 3, 1 3 3, _ 3, e IS.) cos 2 h 1 cos 2. 

+ j ( r ^) cos h l " 2 ( F' > COS 1 2 r 3 1 1 


3(^®Jl)cos h x cos 2 2 1 1 

2 R J 


(A. 3. 4) 


Equation (A. 3. 4) may be used to write down the expression 


for D 


12 ' U 21 


, and D 22 simply by changing the subscripts on 


h and When the double difference t (D u - ° 12 > 

(D - D n is formed, all terms which contain only 
tracking station arguments, K x or h 2 , or which contarn 


t 






only ALSEP arguments, or l 2 , will cancel. Retaining 
only the "mixed" terms of Equation (A. 3 . 4 ) yields 

A4, = |[(Du - D 12 ) - (d 21 - D 22 )] 

= HE { (_5_E) - h 1 ) - cos ( - h 2 ) 

R 

- cos(&2 “ + cos (^2 “ ^2^ ” cos ^1 COS ^1 

+ cos h2 cos + cos cos & 2 " cos h2 cos & 2 ] 

2 

r r 

+ (_!LH) [cos h x cos (£ 1 - h x ) - cos h 2 cos (£ 1 - h 2 ) 

R 

- cos cos(£2 ” h^) + cos h2 cos{&2 “ ^ 2 ^ 

2 

- |cos 2 h 1 cos a l + |cos 2 h 2 cos l 1 + |oos h-j^cos & 2 

2 

IT IT 

- ^cos 2 !^ cos 1 ^] + ( 6 2^ ) [cos cos(&^ - h^) 

R 

- cos a 1 cos (£ 1 - h 2 ) - COS l 2 cos < £ 2 ~ h l ) 

3 2 

+ cos &2 cos (^2 ~ ^2) “ 2 cos cos 

+ ^cos h2 cos 2 + §cos h-^ cos 2 l 2 ~ 2 cos ^2 COS ^ 2 ^ ^ 


(A. 3. 5) 

The expression can be simplified by defining the quantities 
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By substituting in Equation (A. 3. 5) the expressions 
^2 = hp + Ah an< ^ ^2 = ^1 + ^ anc ^ us ^- n 9 the identities 
sin(x ± y) = sin x cos y ± sin y cos x and cos (x ± y) = 
cos x cos y + sin x sin y, we obtain after much algebra the 
result 

4uir r 

0 m 

A<t> = — ^ — [sinA£/2 sinAh/2 cos £ cos h 
1 r e 

+ j(p— )sinA£/2 sinAh(2cos £ cos 2h - sin £ sin 2h) 

1 r m 

+ j(^)sinA£ sinAh/2 (2cos h cos 2£ - sin h sin 2£)] 

(A. 3.6) 


This expression may be rewritten in terms of the parameters 

of Equation (2.2.7) by noting that 2r sin Ah/2 = B , 
r e 

2 (^) sin A£/2 cos £ = Aa, sin Ah = 2sin Ah/2 cos Ah/2 = 

2 

PB/r and cos h = sin (a - A). Making these substitutions 
and making explicit the time dependence of the diurnally vary- 
ing arguments allows Equation (A. 3. 6) to be written as 


Ac}) (t) = a3 ^ Aa { sin ( a —A ( t ) ) + ^[sin 2 (a - A(t)) 


- ^tan £ cos 2 (a - A(t))] 

, r m cos A£/2 cos 2£ r . . 

+ R- S5ST [sin (a - A{t) ) 


- ^-tan 2£ cos (a - A(t))]} 


(A. 3. 7) 
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